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V 

Applications of Perturbative Quantum Chromodynamics 
to Processes with Heavy Quarks 

by 

Alexander Mitov 

Abstract 

In this thesis we apply perturbative QCD to make precision predictions for some 
observables in high-energy processes involving heavy quarks. 

The first apphcation we consider is a prediction for the spectrum of 6-fiavored 
hadrons in top quark decay. For that purpose we calculate at NLO the QCD 
corrections for bottom fragmentation in top decay with the 6-mass fully taken 
into account. Using the perturbative fragmentation function formalism we then 
resum with NLL accuracy large coUinear logs of the ratio of bottom-to-top mass, 
which leads to an essential improvement of the results. Next we perform the 
threshold resummation for the coefficient function for top decay with NLL accu- 
racy. That resummation leads to an important improvement of the b spectrum 
in the region where the produced bottom in top decay carries a large fraction of 
the momentum of the parent top. Finally, we extract information for the non- 
perturbative 6-fragmentation into hadrons from e+e^ data and make a prediction 
for the spectrum of those 6-fiavored hadrons produced in top-quark decay. 



vi 

Our second application is to charm production in charged-current DIS. We 
first calculate with NLL accuracy the soft-gluon resummed coefficient function for 
heavy quark production (initiated by a light quark) in inclusive DIS. Our result is 
applicable for the case of low momentum transfer that is of the order of the mass of 
the heavy quark. We also make a connection of this result to the known result for 
massless quark production. We then apply this result for charm quark production 
at NuTeV and HERA for a wide range of the transferred momentum, and present 
the effect of the threshold resummation on the charm structure functions. 
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Chapter 1 
Introduction 



For the three decades of its existence, the Standard Model (SM) of elementary 
particles has proved to be extremely successful. It has withstood all experimental 
tests and has become a well established theory. All predictions based on the SM 
have been experimentally verified and most of its parameters have been fixed. The 
only sector of the SM that has not been directly experimentally verified as of today 
is the Higgs sector and, in particular, the existence of a neutral massive spin-zero 
particle often simply referred to as Higgs. There is still no direct evidence for its 
existence and, despite the many constraints from precision electro-weak physics, 
the Higgs mass is not known. Not all of the parameters of the Higgs potential arc 
determined, and its Yukawa couplings to the fermions are only implicitly tested 
through the measurements of the masses of the fermions (quarks and leptons). 
Experiments in the near future and most notably the Large Hadron CoUider (LHC) 
and the Tevatron, will either confirm that particle's existence and fix the above 
mentioned parameters, or will significantly increase the limit on its mass. 

The SM describes three of the four presently known interactions: Electromag- 
netic, Weak and Strong. The fourth one - Gravity - will not be considered in this 
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thesis. We will only point out that the orthodox viewpoint - that gravity becomes 
important only at scales of the order of the Plank one and is thus completely ir- 
relevant for present day collider experiments - has been recently challenged. New 
ideas have emerged (Arkani-Hamed, Dimopoulos and Dvali [1] , Randall and Sun- 
drum [2]) suggesting that it may be possible for gravity to play an important role 
at energies as low as a TeV. However, at least at presently accessible energies, the 
experimental data do not favor any gravitational effects (usually through so-called 
Kaluza-Klein modes of gravity or other higher dimensional matter) and most of 
the models deahng with those effects typically assume that the lowest masses of 
the new particles are close to or above 1 TeV. 

The SM is a gauge-field theoretical model based on the following non-abelian 
gauge group: 

GsM = SUi3)c X SUi2)i X C/(l)y. (1.1) 

The gauge group is non-simple and involves three different dimensionless cou- 
pling constants corresponding to each of the three group factors above. The 
SU{2)j X U{1)y part corresponds to the Weinberg-Salam model [3] and provides 
a unified description ^ of the electromagnetic and weak interactions. There are 
four gauge bosons associated with that group: two neutral ones 7 and Z and two 
charged ones W^. The photon 7 is exactly massless due to the unbroken U{1) 
subgroup that is identified with electro-magnetism. The other three are heavy, 
with masses approximately ^ mz — 91 GeV and — 80 GeV; their masses 

arise as a result of electroweak symmetry breaking via the Higgs mechanism. The 
same mechanism also supplies the tree-level masses of the fermions through their 
^Strictly speaking it is not truly unified since the group is not simple and we introduce a 

priori more than one dimensionless coupling. 

^The values of the various constants used throughout this thesis are taken from the Particle 

Data Group [4]. 
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Yukawa couplings to the Higgs doublet. The electroweak vector bosons couple to 
all fermions. The magnitude of those couplings is relatively small; the electro- 
magnetic interactions are suppressed by powers of a — c^/Att ~ 1/137, while the 
effects of the weak interactions are typically proportional to powers of the Fermi 
constant Gp. For that reason the quantum corrections introduced by electro- weak 
interactions are much smaller than the ones due to strong interactions ^. 

The dependence of the dimensionless electroweak gauge couplings g and g' on 
energy is "intuitive" i.e. the strength of the electroweak interactions increase with 
the energy scale. For that reason the SU{2) x U{1) theory of the electroweak 
interactions is a true perturbation theory formulated directly in terms of observ- 
able fields ^. The situation changes dramatically when one considers the strong 
interactions described by the SU{3)c factor in (1.1). First, the coupling constant 
is not small. Moreover, as is well known, a non-abelian gauge theory with such a 
gauge group and with sufficiently small number of active fermions (flavors) exhibits 
the "counter-intuitive" behavior known as asymptotic freedom: the dimensionless 
coupling constant associated with that group decreases with the increasing of the 
energy scale and effectively such a theory behaves as a free theory at high energies. 
However, at small energies the coupling grows and eventually diverges at some fi- 
nite value of the energy scale. That scale, usually denoted as Aqcd, has a value of 
the order of 200 — 300 MeV and quantifies the borderline between the perturbative 
and non-perturbative regimes in such a theory. Although the detailed description 
of QCD will be spared for the next Chapter, from the above discussion it becomes 
clear that such a theory exhibits another remarkable property: confinement. The 

growing of the coupling at low energies (which corresponds to large distances) 

■^This is also the case in the apphcations considered in this thesis as will be detailed in the 
subsequent sections. 

"^With the exception of the quarks; see below. 
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indicates that the particles whose interaction is described by such a gauge theory 
may not be able to exist as free (asymptotic) states at all. Instead, they will form 
bound states. That expectation is confirmed by experiment: no free quarks have 
ever been observed experimentally. Unlike the quarks that carry an additional 
quantum number called color, the observable strongly interacting particles - the 
hadrons - are colorless objects and have the quantum numbers of two-particle 
(mesons) and three-particle (baryons) bound states of quarks. Strictly speaking, 
the property of confinement is an assumption based on the above mentioned be- 
havior of the strong coupling constant and the non-observation of free colored 
particles. The derivation of the properties of the hadrons from the first principles 
(i.e. from QCD) is one of the fundamental problems in theoretical physics today. 

It is a remarkable achievement that we are able to make precise predictions 
for the observed hadronic states in high energy experiments based on a theory 
formulated in terms of non-observable constituents - the quarks. One of the main 
ingredients of the theory that makes this possible is the so-called factorization 
theorem. In essence, it states that in hard scattering experiments with typical 
hard scale Q » h-qcD an observable (e.g. a cross-section) can be written as a 
product (more precisely - convolution) of perturbative and non-perturbative parts. 
The former part can be calculated in perturbation theory based on an expansion in 
the strong coupling constant asiQ"^) while the latter part has to be extracted from 
the experiment. It is possible to extend the factorization theorem even to processes 
where other scales besides Q are present, e.g. when heavy quarks are involved. As 
will be clear from discussions throughout this thesis, there is no absolute notion 
of heavy quarks, i.e. whether a quark is considered heavy or light depends on the 
particular problem being studied. However, in general, light quarks are considered 
to be those with masses below A.qcd (i-e. s) while the ones with masses above 
that scale are usually considered to be heavy (c, 6, t) . The reason behind such a 
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separation is easy to understand: for a heavy quark with mass M, as{M'^) < 1 
and it is indeed sensible to apply perturbation theory. Studies of QCD involving 
processes with heavy quarks are at present an important internal test for QCD 
as well as for obtaining precision predictions that will be needed to distinguish 
signals from new physics (Supersymmetry, Extra Dimensions, etc.). 

QCD has another peculiar feature: there are situations where the convergence 
of the perturbation series is spoiled because of the appearance of additional factors 
that multiply the coupling constant to any order in perturbation series. The 
presence of such terms effectively alters the expansion parameter to a larger value 
which in turn spoils the convergence of the series. To be able to obtain useful 
information in that case, one needs to resum classes of such terms to all orders in 
the coupling constant. 

There are many examples of physical processes where the factorization theo- 
rem plays a decisive role in studies of processes involving heavy quarks and where 
often the application of the above mentioned resummations leads to serious im- 
provement of the perturbative results. This thesis is devoted to the study of two 
such processes in perturbative QCD (pQCD) with detailed phenomenological ap- 
plications: the spectrum of 6-flavored hadrons in top quark decay and neutrino 
production of charm in Deep Inelastic Scattering (DIS). 

In the SM, it is the mass of the top that uniquely distinguishes it from the 
other five flavors. The very large value of this parameter, however, is sufficient 
to draw particular attention to the top quark. Top's large mass is responsible 
for its small lifetime, the latter preventing the top from forming bound states (a 
process known as hadronization) . For that reason, the top quark behaves like a 
real particle and one can safely describe its decay in perturbation theory. Since 
the only experimental information about the top is through its decay products, it 
is very important to have a precise prediction for the decay products of the top. 



CHAPTER 1. INTRODUCTION 



6 



In this thesis we make a prediction for the spectrum of the hadrons resulting from 
the hadronization of the 6-quark in top quark decay and resumming to all orders 
in as and with next-to-leading logarithmic (NLL) accuracy two classes of large 
logs: quasi-coUinear ones that are due to the large ratio of top-to-bottom mass 
and soft ones that arc due to soft-gluon radiation and affect the distribution in 
particular kinematical regions. Such results will be very important after the near 
future high energy experiments supply enough data on top decay. 

The second process studied in this thesis is Deep Inelastic Scattering (DIS), a 
process of scattering of a lepton off a hadronic target. Since the leptons do not 
interact strongly, that process can serve as a tool to study the structure of the 
hadrons by probing them with a virtual vector boson emitted from the lepton. 
DIS is perhaps the most studied high energy process; it was in this process in 
the late 1960's that Bjorken scaling was discovered, and which in turn served as 
a motivation for the introduction of the quark-parton model. Using the available 
next-to-leading (NLO) calculations for heavy quark production in charged-current 
(CC) DIS we calculate the corrections to the charm structure functions in the 
region of low momentum transfer and fully accounting for the mass of the charm 
quark. Such a result is important for the precise determination of the parton 
distribution functions in the target hadrons as well as the value of the charm 
mass. 

This thesis is organized as follows: in the next Chapter we discuss some general 
features of QCD from the perspective of our applications. In Chapter 3 we present 
our original results on 6-fragmentation in top decay. In Chapter 4, we discuss our 
original results for soft-gluon resummation in the coefficient function for heavy 
quark production in CC DIS, which we apply to study the effect of the threshold 
resummation on charm production at low energy transfer. In Chapter 5 are our 
conclusions. In the appendix we have listed some useful results. 



Chapter 2 



Preliminaries: Perturbative QCD 



In the previous Chapter we made some general remarks about the relevance of 
QCD as a theory of the strong interactions as well as few of its peculiar features. 
In this Chapter we are going to systematically review that theory and derive many 
of its properties. The organization of the material in the present Chapter does 
not follow any particular review on QCD but is presented in a way that is suitable 
for our applications. There are many excellent reviews of QCD; some of those 
include [5, 6, 7, 8, 9]. 

2.1 QCD as a Fundamental Model for the Strong 
Interactions 

The strong interactions govern the interactions of hadrons at a wide range of en- 
ergies: from the highest energies accessible to the present day colliders down to 
energies typical for the nuclear physics. At the same time the behavior of the 
strong interaction is very different in the low and high energy regimes. At low 
energies, i.e. energies characterized by a scale « Aqcd, the hadrons behave 

7 
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as if they are fundamental particles. However, no successful field-theoretical de- 
scription in terms of the observed hadrons was found that was able to describe 
the high energy regime and explain the proliferation of observed hadrons at high 
energy colliders. Contrary to the early expectations, the breakthrough towards 
the understanding of the strong interactions was made in the study of the high 
energy behavior of the hadrons. 

2.1.1 Quark Hypothesis 

In 1964 Gell-Mann and Zweig [10] introduced the idea of quarks: a few elemen- 
tary particles that are the building blocks of all hadrons. There are six known 
types of quarks (quark flavors); they are spin 1/2 fermions with rational electric 
charges (in units of the charge of the electron): u.c.t have charge -1-2/3 while 
d,s,b have charge —1/3. The quark hypothesis assumes that the wave function 
of a hadron is constructed from the one-particle wave functions of quarks and/or 
antiquarks. The mesons have the quantum numbers of a quark-antiquark pair 
while the baryons are combinations of three (anti)quarks. Also, in order to avoid 
a problem with the spin-statistics theorem, it was necessary to introduce addi- 
tional hidden quantum number - color [11]. From a comparison with experiment 
it was concluded that each quark flavor must have three different "copies" labelled 
by an additional color index. Since no colored particles have been observed it was 
postulated that the hadrons can only form "colorless" combinations of quarks, i.e. 
color was introduced as an exact global symmetry. It can be shown that the above 
observations plus the requirement that quarks and antiquarks transform under 
different (complex-conjugated) irreducible representations of the color symmetry 
group uniquely fixes the group to be SU{3)c- The quarks transform under the 
fundamental representation 3 of that group, while the antiquarks transform under 
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its conjugated representation 3. 



2.1.2 Parton Model 



The idea that the hadrons are built from elementary constituents - the partons [12] 
- was extremely successful not only in explaining the hadron spectroscopy but also 
in the description of the Bjorken scaling [13] observed in Deep Inelastic Scattering 
(DIS) experiments [14]. The experimental data showed that at large scales the 
structure functions of the nucleons are (approximately) independent of the value of 
the hard energy scale Q and depend only on the Bjorken variable x (for a detailed 
description of DIS we refer the reader to Section (4.1)). The parton model assumes 
that in high energy lepton-nucleon scattering, where the transferred momentum is 
large enough so the masses of the partons and their transverse motion inside the 
nucleon can be neglected, the virtual electroweak vector boson emitted from the 
initial lepton is scattered incoherently by a single free point-like parton a. The 
whole information about the structure of the hadron that is relevant to the high 
energy process is encoded in a scalar function fa{C) called the parton distribution 
function (pdf) representing the probability distribution for finding the parton a 
inside the hadron and carrying a fraction ^ : < ^ < 1 of the momentum of the 
parent hadron. The philosophy of the parton model then suggests the following 
form of the hadronic interactions at high energy: 



where da{h, . . .) is the cross-section for scattering of a hadron h and . . . stands for 
the other particles in the scattering process; da{a . . .) is a parton level cross-section 
with the hadron h replaced by a free parton a, and the partonic cross-section is 
weighted with the distribution function fa{0- Since the parton model is a free 
theory, to lowest order in the electroweak coupling the partonic cross-section is 




a 



(2.1) 
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very simple: da ~ — x), so that the momentum fraction ^ is identified with 
the Bjorken variable x (cf. Section (4.1)). Therefore under the parton model 
assumption (2.1), the structure functions are simply proportional to the pdf fa{x) 
and naturally independent of the hard transferred momentum Q. 

The success of the parton model goes beyond the description of Bjorken scal- 
ing. It also makes a prediction about various relations involving the measurable 
structure functions Fi, F2, . . . (that will be properly introduced in Section (4.1)), 
among which is the Callan-Gross relation [15] : 

F2{x)^xF^{x) , (2.2) 

which is related to the fact that quarks have spin 1/2. The parton model can be 
generalized to other processes as well; one just needs to measure the pdfs for the 
various quarks in a specific process in order to predict a measurable quantity for 

another process. In that procedure, the following relations between the various 
parton distributions following from isospin invariance are often assumed: 

/r"*°"(^) = fr''''"\^) ; fr'"'\^) = fr''"'\x). (2.3) 

2.1.3 QCD: the Dynamical Theory of Color 

Although the parton model was quite successful in the description of many high 
energy processes, it was clear that it is at best a good hint towards a complete 
dynamical theory of the strong interactions. The complete theory would be able 
to explain one of the basic assumptions of the parton model - asymptotic freedom. 
Such a theory was constructed around 1973 after it was understood that the Yang- 
Mills theories do play an important role in high energy physics; at that time the 
renormalizeability of those theories was proved and methods for their quantization 
were developed [16]. It was also shown that the non-abelian theories were the only 
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theories which may exhibit asymptotic freedom, or technically, have negative first 
coefficient in the /9-function (see next section) . All these developments led to the 
construction of QCD as the dynamical theory of the strong interactions [17] as 
follows: 

QCD is a non-abelian gauge theory with six quark flavors. The gauge group can 
be naturally obtained by gauging the exact global color symmetry group SU {S)c- 
The quarks transform under the fundamental representation of SU{3)c- Since 
dim(5'C/(3)) = 8, there are eight gauge bosons called gluons that are electrically 
neutral, carry color charge and, as usual, are hermitian fields that transform under 
the adjoint representation of the gauge group SU{2i)c- 

The lagrangian of QCD has the following form: 

Cqcd = -\f;^F^^^'' + iY.J\l' (D,)^, ^l-T. (2-4) 

where the index i = 1,2,3 runs over the different quark colors and q over the 
quark flavors: q = u,b, s, c, b, t. The field-strengths are given by: 

F;, = d,Al - d^Al - gsr'^AlAl, (2.5) 

and the gauge-covariant derivative is: 

{D,\.^5,,d, + igsY.t1,Al. (2.6) 

a 

In the above equations, gs is the strong coupling constant, /"^"^ are the struc- 
ture constants of the gauge group and are the generators of the fundamental 
representation of the gauge group. In general, the lagrangian (2.4) must be sup- 
plemented with gauge-fixing and ghost terms. The quarks have non-zero masses 
but their origin is outside QCD; in SM their masses result from the electroweak 
symmetry breaking. The only free parameters in QCD are the six quark masses 
and the single gauge coupling constant. 
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The gauge group has the following matrix structure which, for generality, we 
present for arbitrary group SU {N): the fundamental representation has generators 
r, a = 1, . . . , _ 1 that satisfy 

E% = Cf5,,, t,j,k = l,...,N. (2.7) 

a 

For the group SU{3), the generators are usually given by the Gell-Mann ma- 
trices A": = A"/2. Similarly, the adjoint representation has generators T" that 
can be related to the structure constants f"''"^ through: 

(T«);,c = r'", (2.8) 

and 

tr(T'=T'^) = J2 r'^T^'^ = Ca5"^. (2.9) 

a,b 

Above, Cp and Ca are the values of the quadratic Casimir of the gauge algebra 
in the fundamental and adjoint representations respectively: 

- 1 

= Ca^N, (2.10) 

and for the case of SU (3): 

Cf = \] Ca = 3. (2.11) 

Once formulated, it must be shown that QCD indeed is capable of reproducing 
the success of the parton model as a first step. That in fact follows since from the 
very formulation of QCD it is clear that the parton model corresponds to the Born 
approximation of QCD. The real challenge however is to derive the asymptotic 
freedom from first principles and also to derive the corrections to the Bjorken 
scahng. 
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2.1.4 Strong Coupling Constant 

The running of the renormahzed strong couphng as = Qs/'^'^ is determined from 
the following equation: 

= 2P{as). (2.12) 

The /3-function I3{as) ^ has a series expansion in the coupling constant as- It can 
be determined up to a fixed order in perturbation theory from explicit evaluation 
of the gauge coupling renormalization constant Zg-. 

At present, the /5-function of QCD is known to four loops in the MS scheme [18]. 
However, since for all applications in this Thesis we need the evolution of the 
strong coupling to two loops, we will present only the two loop result: 

I3{as) = -boal - ha^ (2.14) 

with bo and bi given by 

33 -2n/ , 153-19nf ^ 

The first two coefficients of the /3-function are independent of the renormalization 
scheme. However, that is no longer true for the higher order terms. In Eq.(2.15), 
71/ is the number of active massless flavors. In the presence of quark masses, its 
value becomes scale dependent. If one assumes strong ordering among the quark 
masses, then for scale /i : m„ <<//<< m^+i all flavors with masses below m„+i 
are effectively massless while the rest of the flavors are heavy and can be integrated 
out. An example is the case of top decay considered in Chapter 3 where the scale 
^In contrast to most of the standard presentations (e.g. [5]), we have introduced P{as) 
through the relation: (3{as) = {gs/^'i^)(3{g). 
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is running between the b and the t masses and because of their large difference 
we simply fix n/ = 5. However, in practical applications (especially with scales of 
the order of the b and the c quark) the strong ordering assumption is not always 
valid and therefore the choice of n/ is somehow ambiguous. The common practice 
is to change the value of ^/ by one unit when the hard scale crosses the mass of 
the corresponding heavy quark. Such a change should be supplemented with an 
additional constraint that relates the values of the strong coupling evaluated in 
the two schemes at the switching point. In the MS renormalization scheme, the 
strong coupling is continuous at the switching points [19] ( see also [4]), up to 
negligible corrections of order 0{ag). 

Now it is easy to show that indeed QCD enjoys the property of asymptotic 
freedom [20]. In a regime where the strong coupling is small, from (2.15) and 
(2.12), it is easy to see that the strong couphng is a decreasing function of the 
scale /X if the number of flavors n/ < 33/2. That requirement is satisfled in QCD. 
The exact solution of Eq.(2.12) to NLO is given by: 



I [ 6iln[ln(/xVA ^)] 

6oln(//VA2) [ bl\TL{iJ?/K^ 



One can use this expression in order to relate the values of the strong coupling at 
two different scales with NLO accuracy [21]: 

' l + b^as{^l^)Hky|l^)\ boi + boas{pi')Hkyi,^) ' 

X ln(l + boasii^') HkVpi')) + O (4(/i2)[a5(/x^) H^V i^'T) ) ■ 

The constant A contains all the information about the boundary condition to 
which Eq.(2.12) must be subjected. It is a low energy scale where the strong 
coupling diverges. As we mentioned in Chapter 1, A represents the border between 
the perturbative and non-perturbative regimes of QCD. In practice the value of A 
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is not unambiguous. In high energy experiments one typically obtains information 
about the strong coupling constant at some large scale and only from there the 
value of A is inferred. It is obvious however, that in this way the determination of 
A "absorbs" all ambiguities such as the dependence on the order of as (LO, NLO, 
etc.), the choice of the value of and the choice of the renormalization scheme. 
In this Thesis we use the following value of the strong coupling at NLO [4] : 

a^(M|)= 0.118 (2.18) 

It leads to A^^) ^ 200 MeV and A^^) ^ 300 MeV. The precise values used in our 
applications will be discussed in the next Chapters. 

2.2 QCD Factorization Theorem 

QCD is formulated in terms of quarks and gluons while the experimentally ob- 
served states are hadrons. Since at present we are not able to describe the non- 
perturbative regime of QCD, a fruitful application of QCD to the hadronic in- 
teractions requires a universal way of splitting the contributions from short- and 
long-distance physics. 

Such a separation is possible. It is known as the (QCD) factorization theorem 
[22, 23] and states that for processes that have initial and/or observed final state 
hadrons the differential cross-section has the following form: 

h,h' ij 

(g) Dh'/f{x,fj,'') + 0{A/Q) (2.19) 

The factor fi/h stands for the parton distribution function of the parton i inside 
the hadron h present in the initial state, and Q and x represent the hard scale and 
some kinematical variable respectively. Unlike the simple parton model (2.1), the 
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parton distributions also depend on the factorization/renormalization scale /x^. 
That difference has important consequences since it incorporates the violations of 
Bjorken scaling. The second factor dui-^f, also known as the (Wilson) coefficient 
function, represents the partonic hard scattering cross section for the reaction 
i ^ f that depends on the unphysical renormalization and factorization scales /i^ 
and and on the masses of the heavy quarks m^. The last factor in (2.19) is the so 
called fragmentation function D. It contains the information for the hadronization 
of the hard parton / (that is produced in the hard process described by the 
partonic cross-section da) into an observed hadron h' . The integral convolution 
appearing in (2.19) is defined as: 

/■I rl 2 

U®9){x)= / -f{z)g{x/z), (2.20) 

Jx Z 

where / and g are two functions with argument x : < x < 1. The convolution 
has the property (2.40) which is important in practical applications. 

The real power of the factorization theorem is in the fact that the distri- 
bution/fragmentation functions are universal: they depend only on the non- 
perturbative transition they describe and not on the hard scattering process. 
That is why once they are measured in one process, they can be applied to any 
other process. At the same time, the coefficient function contains all the infor- 
mation about the hard scattering process and is independent of the details of the 
non-perturbative transitions. Although the (8)-product of coefficient function and 
distribution/fragmentation function is an observable and therefore free from any 
ambiguity, the distribution, fragmentation and the coefficient functions are sep- 
arately ambiguous. In particular, they are scheme dependent; the origin of that 
scheme dependence is in the treatment of the IR divergences associated with their 
computation. Let us describe that in more detail: 

The evaluation of the coefficient function proceeds in the following way: one 
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calculates the hadronic process da that is under study by formally replacing each 
initial hadron h with an on-shell parton a. To that end, and in accordance with 
the factorization theorem (2.19) we just described, one introduces (also formally) 
new parton distributions fi/a which have the meaning of a distribution of a parton 
i inside a parton a; wc treat the fragmenting partons in a similar fashion. The 
main purpose of the fictitious distributions fi/a is to absorb all the IR sensitive 
contributions from the calculated cross-section. As a next step, one simply dis- 
cards the functions fi/a and what is left is the needed coefficient function da. The 
extraction of the partonic pdfs is physically equivalent to the absorbtion of the 
IR sensitive contributions into the pdfs fi/h (see also Section (2.4.1)). 

Clearly, such a procedure is very similar (and is in fact related) to the UV 
renormalization where one introduces appropriate counter-terms to absorb (and 
thus cancel) the UV divergences appearing in the Feynman diagrams. In the most 
often used MS subtraction scheme the partonic pdfs read: 

fi/a{x) = M(l - ^) + ^ (-7 + 7^ - 1^(4^)) + 0{al), (2.21) 

where Pij\x) are the leading order AltareUi-Parisi splitting functions that will 
be defined and thoroughly discussed in the next section, and e — (4 — D)/2 (see 
the appendix for more details). The subtraction scheme for the IR divergences 
emerging in the limit e ^ is related to the renormalization scheme used to 
remove the UV divergences appearing in the formal (operator) definition of the 
parton densities. Following the approach in [22], the quark distribution can be 
defined (in the hght-cone gauge = 0) as: 

Mx) = ^J dy-e-^^^^^-(/i(P)|V^(0,y-,0^)7+V'(0,0,0^)|/i(P)), (2.22) 

with similar results for the antiquark and gluon distributions. In Eq.(2.22) P is 
the momentum of the hadronic state h, — 1/V^{V^ ± V^) is the light-cone 
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co-ordinate of a four- vector V and ijj is the operator of the quark field that carries 
momentum fraction x : — xP^ . For a partonic initial state, one can evaluate 
the partonic PDF in powers of as by explicit calculation of the matrix element in 
(2.22). The result to order as shows the following structure of the singularities: 

U,a{x, e) = 5,„(5(1 - x) + ^ (- -) pII\x) + UV counterterm. (2.23) 

iTT \euv ^iRy 

Since by the UV renormalization the UV counterterm cancels the 1 / euv pole, we 
have a remaining IR divergence in the partonic PDF. The explicit form of that 

singularity obviously depends on the chosen UV renormalization scheme; for the 
MS scheme the result is presented in (2.21). Similar considerations apply for the 
distribution functions as well; that is discussed e.g. in [24] and [25]. 

The factorization theorem in the presence of massive quarks is accurate up to 
terms 0(A/Q); see Eq.(2.19). The remainder in (2.19) contains terms that are 
power suppressed and are uniform in Q i.e. no terms with suppression ~ {m/Q)"' 
are present [23]. The proof of the factorization theorem presented in [23] uses a 
variable flavor number scheme (VFNS). The VFNS treats the light quarks u, d and 
s as massless and always includes them as active flavors in the running of the strong 
coupling as- The treatment of the heavy flavors c, b and t is, however, process 
dependent. If the typical energy scale is below the corresponding quark mass 
then that quark is treated as heavy and is integrated out. In particular it does not 
contribute to the evolution of the strong coupling and does not have an associated 
parton density. The quarks with masses below the hard scale are treated in 
a different way: they contribute to the strong coupling as if they are exactly 
massless, and they have their own distribution functions ^ which are evolved with 
^ Those arc introduced in order to systematically resum large logs of collinear origin that 
appear to all orders in as- Schemes without heavy quark densities exist and are called Fixed 
Flavor Number Schemes (FFNS). An example is the GRV 98 set of parton distributions [26]. 
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the energy scale via evolution equations with massless kernels (see next section). 
Such a scheme corresponds to the renormalization scheme of CoUins, Wilczek 
and Zee [27] and is also used (through the so-called Aivazis-CoUins-Olness-Tung 
(ACOT) prescription [28]) in the evolution of the CTEQ parton densities [29]. 

The physical picture behind the factorization theorem is quite simple. One for- 
mally introduces a scale /ip which separates the short- from long-distance physics 
involved in the process. It is intuitively clear that such separation must indeed 
occur in the hmit of large values of the hard scale Q. The time scale for the hard 
interaction is of the order of and therefore quite small, while the typical time 
for the hadronization effects is not smaller than Aqqjj ^ Q^^. As a result, in the 
limit (5 — > oo, short- and long-distance effects cannot interfere with each other 
and therefore factorize. More formally, the separation between small and large 
scales means that all contributing Feynman diagrams that have lines with small 
virtuality can be separated from the lines with large (of the order of the hard 
scale Q) virtuality. The former diagrams constitute the distribution functions 
while the latter give the hard coefficient function. Such non-trivial factorization 
for the terms with leading power in l/Q^ (the so-called leading twist terms) was 
proved by Libby and Sterman [30]. If one works in a covariant gauge then there 
are present additional gluon hues that connect the small- and large-scale factors. 
However, those additional contributions cancel when summed up in a gauge in- 
variant subset of diagrams. These contributions disappear, however, if one works 
in a physical gauge n.A — 0, where n is some fixed four- vector. The physical 
gauge also has the important property that each diagram that contributes to the 
physical cross-section can be given a direct partonic interpretation [8] . 
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2.3 Perturbative Evolution: DGLAP Equations 

In this section we will devote our attention to the dependence of the various factors 
in Eq.(2.19) on the renormalization and factorization scales which incorporate the 
scaling- violation effects. To better illuminate our point, we are going to make the 
following two simplifications throughout this section: first, we will set /ir — /J'F = 
jji. This is a standard choice in the studies of pQCD which, however, will not 
restrict the generality of our discussion. If needed, the separate dependence on 
both scales can be easily restored with the use of the running of the strong coupling 
(see Eq.(2.17)). The second simplification is that we will consider Eq.(2.19) as 
having a single fragmentation or distribution function multiplying the coefficient 
function da. We will consider those two "representative" cases (only initial or 
final observed hadrons) of Eq.(2.19) separately. 
We start with: 

da{x,Q'^,m^) =Y, fi/N{x, /i'^) dai^x{x,Q^,m^, n^), (2.24) 

i 

which corresponds to the case of a single hadron (nucleon N) in the initial state 
and no observed hadrons in the final state. A prominent example is the case of 
inclusive Neutral Current (NC) or Charged Current (CC) DIS: 

e + N ^ £' + X, (2.25) 

with i and i' being leptons, a hadron (usually a nucleon) and X stands for any 
unobserved hadrons produced in the reaction (2.25). This case describes reactions 
with so-called space-like evolution. 

As a representative for a reaction with a single fragmentation function we take: 

da{x,Q'^,m^) dae+e--^f{x,Q^,m?, f/) (g) Dh/f{x,f/), (2.26) 

/ 
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which corresponds to the case of inclusive production of a single hadron /i in a 
non-hadronic collision. These reactions are known as having time-like evolution. 
An example is the case of a single particle inclusive e'^e~ annihilation: 



with h being an observed hadron and, as usual, X stands for any unobserved 
hadrons produced in the reaction (2.27). Also, for brevity, we have omitted the 
remainders in Eqns.(2.24) and (2.26). 

2.3.1 The Case of Space-like Evolution 

Let us concentrate on Eq.(2.24). Since the LHS is independent of /i we can set 
that scale to any value we like. Among all the possible values, the choice /x^ = Q'^ 
is particularly appealing as will become clear below. 

From dimensional considerations the coefficient function can be written as 



where for the time being we consider the case when no heavy quarks are present; we 
will generalize our considerations in the next section. In Eq.(2.28) as is the Born 
cross-section for the partonic subprocess and the function C is a dimensionless 
function that has a power series decomposition in the strong coupling as- 

It is now obvious that upon setting /i^ — the coefficient function takes 
the form C{1, as{Q^)) and depends only on the strong coupling and on no other 
large (or small) parameters. Since as is evaluated at the large scale where 
the former is small, the coefficient function can be easily and efficiently calculated 
to some fixed order in perturbation theory. However, as a result of the choice of 
scale we have made, the distribution function has now become Q-dependent. That 
dependence is very important. It indicates that the universality of a distribution 



e"^ + e 



h + X 



(2.27) 




(2.28) 
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function may be reduced since the PDF is specific to the experimental setup where 
it is extracted and therefore cannot be apphed to processes with different hard 
scale. 

Fortunately, there exists a way to relate distribution functions at different 
scales. The scale dependence of the fragmentation functions /j, i = q,q, g is pcr- 
turbatively controlled and is given as a solution to a system of integro-differential 
equations known as Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equa- 
tions [31, 32]. That way the universality of the distribution function is retained; 
we only need to extract from experiment the distribution functions at one given 
scale Qq. Then that input can be used as the initial condition for the DGLAP 
equations and the PDF at any other scale can be predicted. In practice that pro- 
cedure works in the following way: at some low scale Qo ~ 1 GeV one writes down 
a function of z, that contains small number of free parameters. Then one evolves 
that initial condition via the DGLAP equations to different scales where experi- 
mental data exist, and one tries to fit those data by adjusting the parameters of 
the initial condition. 

The DGLAP equations are: 



and describe in general a system of 2n/ -|- 1 equations for the distribution func- 
tions of all flavors of quarks, antiquarks and the gluon. The kernels have 
perturbative expansions in powers of the strong coupling: 



Pij (z) are the Altarelli-Parisi splitting functions [31] that also appeared in Eq.(2.21). 
Because color and flavor commute, those functions are independent of the quark 




(2.29) 




(2.30) 
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flavor i.e. P^^^.{z) = 5ijP^{z). They also satisfy other relations as a result of 
probability conservation: 

£dzPi^\z)^0, (2.31) 

and momentum conservation: 

£dzz{P^^\z) + P^^\z))=0, 

£ dz z (2n,p(J)(z) + P(;)(z)) = 0. (2.32) 

Note that because of the property (2.31), the functions P^q\z) are not positive 
definite. They are distributions instead. The kernels Pij satisfy also the follow- 
ing important relations as a result of charge invariance and the SU{nf) flavor 
symmetry: 

P = P- - ■ P - = P- 

^ng ~ ~ ^gg ' ^gn ~ ^ggi ~ ^gg- (2.33) 

The explicit expressions for the four independent leading order splitting functions 
are: 



+ H^l^id-.). (2.34) 

The distribution (F(^))+ is defined as S^{F{z))^H{z)dz = /q F{z){H{z)-H{l))dz, 
and its properties are discussed in the appendix. 
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One can substantially simplify the study of the DGLAP equations if one takes 
advantage of the flavor symmetry. As we mentioned earlier, the DGLAP equations 
describe the evolution of 2n/ + 1 partons that are all massless (see also the discus- 
sion following Eq.(2.23)). Prom the form of the QCD lagrangian with Uf massless 
quarks (2.4), it is evident that the theory has an additional global SU{nf ) flavor 
symmetry. Since the gluons are flavor-neutral, they transform as singlets under 
the flavor group. The quarks in general transform nontrivially under that group. 
One can split the 2nf (anti)quark fields into one singlet: 



and 2n/ — 1 non-singlet (NS) combinations, n/ of the NS fields can be taken as 
the differences M~ — fqi{x,ii) — fq.{x,iJ,). The other rif — 1 combinations, which 
we denote by Mj~, depend on the value of n/ and can be found in [6]. Clearly, the 
NS combinations do not mix with the singlet; in particular they do not mix with 
the gluon. To LO, all NS fields also spht from each other so we have a separate 
equation for each NS field: 

with P^^^{z) given in Eq.(2.34). Beyond the leading order, however, the evolution 
kernels are no longer flavor diagonal. One can still write the evolution equations in 
diagonal form that is similar to the LO case Eq.(2.35), but the kernels P^P' and 
^NS corresponding to the fields and M~ are now different. Their explicit 
expressions can be found in [6] as well. 

In the singlet sector, there is non-trivial mixing between the gluon density 
g and the quark singlet state E. The kernels of the evolution equations for the 
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"two-vector" {Hlz, /j,), g{z, /j,)) form a 2 x 2 matrix: 

^ Pqq{z,asin'^)) 2nfPqg{z,as{n'^)) ^ 

\ Pgg{z,as{fi'^)) Pgg{z,as{fi^)) J 

The NLO kernels in the singlet sector can be found in [6]. Similarly to the NS case, 
we do not present them here because of their length. The original derivations are 
presented in [25, 33]. For future reference we will only present the large z behavior 
of the MS sphtting functions at NLO: 

lim P„ = C.f (l + + 0(a|)) ^ (2.36) 

where: 

The origin of the DGLAP equations is in rcnormalization group invariance. 
That invariance is manifested as independence of physical quantities (for example 
a cross-section) of the renormalization scale /ir that is introduced as a result of 
the renormalization procedure. Eq.(2.24) is a typical example. Before we exploit 
the fact that the LHS is independent of fx, let us take the MeUin transform of 
Eq.(2.24) (with all masses omitted for now): 

daiN, Q') = Yl fii^^ l^^)dai^xiN, Q^ /z^), (2.38) 

i 

where the Mellin transform F{N) of a function F{z) is defined as 

F{N) = f^dz z^-^F(z). (2.39) 
Jo 

To derive Eq.(2.38) we used the following important property: the Mellin trans- 
form of a convolution of two functions is the ordinary product of the Mellin trans- 
forms of those functions, i.e. 

dz z^-\F G){z) = F{N)G{N). (2.40) 



/ 

Jo 
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Taking /^-derivatives of both sides of Eq.(2.38) and omitting the summation 
over the partons (as in the NS case for example), one has: 

/.|^ln/(7V,/.2) = -l^-^lnC (n, ^,c.s{i^')^ = 2^{N,as{i^')). (2.41) 

In Eq.(2.41) we have also used Eq.(2.28). From the derivation of Eq.(2.41) it 
is clear that the function ^{N,as{fJ''^)) plays the role of a "separation constant" 
that depends on the common variables of the functions / and C. Therefore, that 
equation is of importance only if the function 7 is known. Let us assume for the 
moment that this is the case; then with simple integration one can obtain the 
scale dependence of the moments of the distribution function /: 

m /.^) = m exp J^^ ^j{N, as{e))j , (2.42) 

if the value of the fragmentation function at some initial scale /iq is known. Com- 
bining this result with Eq.(2.24) and after setting fi — Q we can finally complete 
the factorization program discussed in the previous section. To calculate a phys- 
ical observable we only need the PDF at some initial scale and the coefficient 
function for the hard process which is computable in perturbation theory since it 
does not contain any large parameters: 

da{N,Q) = aBf{N,,,l)C{N,l,as{Q')) exp H —j{N,as{P))j . (2.43) 

Of course, in practical applications one needs to carefully account for the mixing 
between the various distribution functions that would appear in the relation (2.43). 

To make contact with the DGLAP equations, we need only notice that by 
setting 

7.,(iV, asifi')) = I' dz z^'-'P^.iz, asifi')), (2.44) 

and applying a Mellin transformation to both sides of the DGLAP equation we 
get exactly the equation (2.41) for the fragmentation function /. That way we 
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conclude that the functions 7jj are the Melhn moments of the evolution kernels 
P 

The functions indeed have a deep physical meaning which is not hard to 
understand. For that purpose let us rewrite schematically Eq.(2.41) for /: 

+ 27(^5)) In/ = (2.45) 

Eq.(2.45) resembles the renormalization-group equation of a matrix element of 
an operator in a massless theory with anomalous dimension 7. The functions 
•yiji^N, as{^i'^)) are indeed anomalous dimensions but of local operators of spin-N 
and twist-two that appear in the light-cone operator product expansion (OPE) of 
two currents. The derivation of this result is rather lengthy and can be found in 
[34, 35]; for a more pedagogical discussion see [5]. Here we shall only mark the 
main steps in the example of DIS. The DIS cross-section is of the form L^^^W^i, 
where L is a tensor constructed from the lepton's four-momentum and W has the 
form: 

W,, ~ j d''ye^^-y{P\[J^{y), J.(0)]|P). (2.46) 

The non-local product of the two electromagnetic currents J^{y)J^{0) is sand- 
wiched between the hadronic state \h{P)) and its conjugate. This operator prod- 
uct can be expanded in a series over local operators 0„: 

J{y)JW = Ec'niy'mr-vr^, . . (2.47) 

n,i 

Sandwiching the current product in the hadronic states: 

{P\0!;,y"\P) = AnA^^')P'' ■■■P"", (2.48) 
results in the following expression: 

F(n, Q') = A„,,(/.2)C;(g7/.^ as). (2.49) 
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Above, F is the n-th Mellin moment of any of the structure functions. Eq.(2.47) 
imphes that the coefficient function C satisfies the renormahzation group equation 
with anomalous dimension 277—70- Since the current J is conserved its anomalous 
dimension vanishes. The anomalous dimension corresponding to an operator O 
can be computed in perturbation theory from the relation: 

7o = f|lnXo. (2.50) 

where Zq is the renormahzation constant for the operator O. 

Collecting the results above, we see that after setting fj, — Q we obtain 
Eq.(2.43) for the DIS structure functions, and the anomalous dimension 7(A'") 
presented there corresponds to the local operators Oj^. 

Let us make few more brief comments on the DGLAP equations. First, the 
DGLAP equations in the form that we discussed were first introduced by Altarelli 
and Parisi in [31]. There the authors were following a more or less reversed 
approach: starting from the OPE results we discussed above, they introduced the 
AP sphtting functions through the relations (2.44). The important result is that 
not only arc the 2;-spacc evolution equations equivalent to the renormahzation 
group ones, but the former offer new way of deriving the splitting functions based 
on their partonic interpretation. Second, similarly to the beta-function, the one- 
loop splitting functions Pl^j\z) are renormahzation scheme independent. However, 
that is not true for the higher order functions which are renormahzation scheme 
dependent. The most common choice is to work in MS scheme. 

2.3.2 The Case of Time-like Evolution 

All the considerations that were made for the case of Eq.(2.24) can also be made 
for the fragmentation case (2.26). There are, however, a few differences between 
those two cases and we will discuss them now. 
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The fragmentation functions D have an interpretation different from that of 
the parton distribution functions. The function Di/h{z,ii'^) represents the proba- 
bihty density that a parton i produced at scale //^ will fragment to an observed 
(and therefore on-mass-shell) hadron h. Similarly to the distribution functions, 
the evolution of the fragmentation functions is also described by the DGLAP 
equations. The one- loop splitting functions Pfime-Hke ijl-^) coincide with those in 
the space-hke case (2.34). However the time-hke and the space-hke evolution ker- 
nels differ beyond the leading order. The NLO time-hke functions can be found in 
[25, 33]. The large z behavior of the NLO time-hke evolution kernels is the same 
as for the space-like kernels, Eq.(2.36). 

2.4 Infrared Effects 

In our previous discussion we neglected the presence of masses of the quarks. 
However for a theory which is sensitive to the IR such as QCD, a detailed account 
for those effects is needed. 

It is well known that in a gauge theory with massless matter fields, in addition 
to the usual UV divergences, there is another type of divergences that occurs 
in the evaluation of the Feynman diagrams. These are known as infrared (IR) 
divergences and as shown by Sterman for QCD in [36] can be divided into two 
types: coUinear and soft. 

The coUinear divergences are due to the vanishing mass of the radiating par- 
ticle (usually the quarks) . When a quark radiates a gluon that is almost coUinear 
to it, the corresponding real or virtual emission diagrams diverge. One can reg- 
ulate such a divergence by introducing a small quark mass or by working in 
D = (4 — e)/2 space-time dimensions. Then a coUinear divergence shows up 
as logarithmic singularity ~ In(m^) or as 1/e pole respectively. In principle, since 
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quarks have non-vanishing masses, the quantities calculated should be finite and 
free of IR coUinear singularities. That is not the case, however. In the pertur- 
bative regime it is not the absolute value of the quark that is important but its 
value with respect to some typical scale. If that typical scale - usually the hard 
scale Q - is much larger than the quark mass, then in a pcrturbative calculation 
there appear large logs ln{m^/Q^). Although these logs are finite, they appear 
to any order in perturbation theory and systematically multiply the strong cou- 
pling constant. Thus, the effective perturbation parameter is not as anymore but 
as X (a power of Inm^/Q^). The latter can be quite large and can even invalidate 
the perturbation series. In effect, small but non-zero quark mass leaves the result 
finite but unphysical; one should sum up to all orders terms of this type in order 
to be able to make definite perturbative predictions. Such large logs are called 
quasi-coUinear logs and are classified in the following way: a term at order has 
the form 



Terms with k — n are known as leading logarithmic (LL) terms, the ones with 
k — n — 1 are the next-to-leading logarithms (NLL) etc. It is a peculiar feature 
of QCD that due to its non-abelian gauge group not only quarks but also gluons 
can radiate coUinear gluons. Unlike the quarks however the gluons are exactly 
massless due to the gauge symmetry. 

The origin of the soft divergences is in the vanishing mass of the gauge fields 
(the gluon). These divergences manifest themselves as singularities in the loop 
integrals over gauge boson lines in the kinematical region where the energy Uq of 
the gluon vanishes or in the real emission diagrams where gluons with vanishing 
energy are emitted. The most convenient way to regularize those divergences is to 
work with dimensional regularization since it preserves both Poincare and gauge 




(2.51) 
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invariance (unlike gluon-mass regularization) . 

It was understood long ago [37] that the problem of IR divergences is rooted 
in the way the physical observables are defined. It is intuitively clear that a 
state containing a hard parton cannot be distinguished from a state containing in 
addition arbitrary number of soft (or collinear) gluons. In those singular limits 
the particle nature of the real soft (or collinear) gluon is not well defined and as 
a result we need to deal with degenerate states. The conditions for cancellation 
of the IR divergences are stated in the following theorem: 

Kinoshita-Lee-Nauenberg Theorem [38]: In a theory with massless fields, 
transition rates are free of IR divergences if a summation over the initial and final 
degenerate states is carried out. 

The proof of the theorem can be found in [5]. Its content is however clear: 
a physical state is one that contains arbitrary number of soft (collinear) gluons. 
When applied to pcrturbative calculations, the KLN theorem means that to some 
fixed order in the coupling constant one should take into account the contributions 
from virtual and real emission diagrams with arbitrary numbers of radiated gluons. 
Only their sum will be IR finite. That way we arrive at the idea of an inclusive 
observable: a calculated cross-section will be IR finite if it does not distinguish a 
state with one particle from a state with a number of soft (collinear) gluons. We 
will see examples in the next Chapter when wc discuss the decay of the top quark. 

In terms of Feynman diagrams (to all orders) , the IR divergences are generated 
only from real or virtual emission lines connected exclusively to external (hard) 
lines in the diagram. The reason is that the internal lines are typically highly off- 
shell and thus regulate any possible divergence. That observation is important for 
constructing an explicit proof of cancellation of the soft divergences and was used 
first by Weinberg [39] in the context of QED. Let us also mention that in fact that 
property leads in the context of QCD to the factorization of the IR singularities 
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in the hard diagrams. This is a very important property that we will make use of 
in the next section. For an excellent discussion see [40] . 

2.4.1 Heavy Quark Masses 

Let us return to Eq.(2.24) or (2.26) and now take into account the masses of the 
quarks that we neglected in the discussion in the previous section. For example, 
Eq.(2.28) now generalizes to: 



It is clear that whatever choice we make for the scale /x, we cannot set equal to 
one both mass ratios that appear in the RHS of Eq.(2.52). Just for the sake of 
being less abstract, let us set as before /s^ — Q^. Then (2.52) takes the form: 



Since the ratio m^/ can take any non- negative real value, we anticipate a strong 
dependence on the value of the quark mass. The case when ^ was already 
discussed in Section (2.2): one can simply integrate out the heavy quark and work 
in an effective theory where that flavor is omitted. As a result, the deflnition of 
the coupling constant becomes dependent on the number of flavors lighter than 
Q and satisfies non-trivial conditions at the switching points. To understand the 
dependence on in the case when is not (much) larger than Q^, we first need 
to know if the cross-section da is IR safe, i.e. if it is finite in the hmit — > 0. If 
it is coUinearly safe, we can represent it as: 



Such cross-sections are well behaved and can be obtained by explicit calculation 
to any order in perturbation theory. A typical example is the case when da is a 




(2.52) 




(2.53) 




(2.54) 



total partonic cross-section, e.g. for the process e+e 



hadrons. 
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The case when da is not IR safe is more comphcated and at the same time 
perhaps more common. In this case the hmit m — > is singular, i.e. da di- 
verges. Examples are the cases where da is an inclusive differential cross-section 
for production of a parton, or a process which is initiated by a single parton 
(the generalization to multiple partons is straightforward). According to (2.19), 
the hadron level result is a convolution of the partonic cross-section with distribu- 
tion/fragmentation function. However in the presence of IR divergences we should 
first understand how to make sense of such divergent results. 

Let us first consider the case when the parton in consideration is light i.e. we 
take ~ 0. The first thing to note is that since it is experimentally measurable, 
the physical process of creation of a hadron or the process that is initiated by a sin- 
gle hadron is not IR divergent. Such a process undergoes complicated stages that 
are controlled by long-distance physics. However, our present ignorance about 
confinement makes us simplify the calculation by assuming that the produced 
hadron has been created from a single parton, which later non-perturbatively 
hadronizes (the latter process being described by the fragmentation function D). 
Such a description is typical for production of light hadrons (i.e. hadrons that 
are constructed out of light quarks) . Thus in simplifying the process we necessar- 
ily introduce mass singularities. The understanding of their origin suggests the 
method to cure them: one calculates the partonic cross-section in perturbation 
theory that is (usually) regulated in dimensional regularization. As we discussed 
in Section (2.2) those divergences factorize and can be subtracted in a particular 
scheme (usually MS). Then the subtracted partonic cross-section is convoluted 
with a non-perturbative fragmentation function which is process independent but 
subtraction scheme dependent. The same considerations apply for a process with 
space-like evolution. The justification for such a procedure is that the subtraction 
is physically equivalent to the absorbtion of the effects sensitive to long-distance 
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physics into the distribution/fragmentation functions. 

After the subtraction of the mass singularities and defining the correspond- 
ing distribution/fragmentation functions we are in the situation described by 
Eq.(2.43): 

da{N,Q)^aBC{N,l,as{Q')) expU^^ —^{N,as{k'))jD{N,i,l) (2.55) 

where we again consider the non-singlet case for simplicity. The function C here 
describes the (MS) subtracted partonic cross-section while D is an initial condi- 
tion for the fragmentation function that is extracted from experiment where light 
hadrons are measured at some low scale /iq- 

2.4.2 Perturbative Fragmentation Function Formalism 

Let us turn our attention to the case of coUinearly-sensitive processes that involve 
heavy quarks and typical hard scale Q somewhat larger than the quark mass. As 
usual, by heavy we mean c,b or t quarks. Although the partonic cross-section 
for such processes is divergent in the zero- mass limit, we cannot really set the 
masses of the heavy quarks to zero since they are not that small. We are then 
in a situation that we previously described: the results are formally finite but in 
practice perturbation theory cannot be applied in a straightforward way since, as 
shown in Eq.(2.51), potentially large logs ln(m^/(5^) appear to all orders in as- 
There are two physically different cases where such a situation can occur. 

The first is the case of partonic processes that are initiated by quarks with 
non-zero mass. Although this is an interesting problem, it is outside the scope of 
this Thesis since it mostly concerns the non-trivial treatment of heavy flavors in 
the nucleons. We will only indicate that in such situations one typically would 
subtract the divergent part of the coefficient function and then convolute it with 
the usual massless parton density. An example is the case of [41] where the strange 
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quark was treated as having non-zero mass and after the subtraction of the quasi- 
coUinear logs the effect of the (finite) power corrections of the mass of the strange 
quark was studied. 

The second case is the one where a heavy quark (usually c or b) is created in 
a hard collision with a typical scale Q. Such processes contribute to the creation 
of c- or 6-flavored hadrons and can be described within the formalism of the per- 
turbative fragmentation function (PFF) [42]. This approach has been extensively 
used for e'^e~ annihilation [43, 44, 45, 46, 47], hadron coUisions [48, 49] and pho- 
toproduction [50, 45]. In the next Chapter we describe its application for bottom 
quark production in top quark decay t bW [51]. 

According to the factorization theorem (2.19), we can write the cross-section 
for creating a heavy-flavored hadron H in the following way: 



where (cr*) is the cross-section for production of hadron H (heavy quark q) 
and Dnp is a non-perturbative function that describes the transition g — > if at 
the scale set by the mass m of the heavy quark q. That function is to be obtained 
from a comparison with experiment. The kinematical variable z describes the pa- 
rameter (s) of the observed final state and typically is an appropriately normalized 
energy fraction. The normalization of the cross-sections is chosen in such a way 
that the Born term equals exactly 5{1 — z). 

The partonic cross-section cr* can be calculated in perturbation theory. It is 
finite because of the finite mass m but is not IR safe since the process is not 
completely inclusive. It contains large quasi-collinear logs to all orders in 0:5 that 
must be resummed as we previously discussed. The resummation of those logs is 




(2.56) 
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done in complete analogy with Eq.(2.55). One writes: 



1 dai 
(Jo dz 




(2.57) 



where a is the cross-section for production of massless parton i with the coUinear 
singularity subtracted in the MS scheme. The function Di{z, 11,171) is the PFF 
and it describes the transition of a massless parton i to a massive quark q. 

The ansatz (2.57) has the following physical interpretation: in a hard collision 
set by the large scale Q a massive parton is produced at large transverse momen- 
tum. For that reason it behaves like a massless parton. The replacement of a 
massive parton with a massless one after the collinear divergence is subtracted is 
justified up to powers of rn^/Q^. The formalism of the PFF is applicable when 
such power corrections are small and can be neglected. For that reason, the cross- 
section a is insensitive to low energy (i.e. at the order of m) physics and depends 
only on Q but not on m. The scale fi is the factorization scale that, as usual, sep- 
arates the low from the high energy regimes. Similarly, the function Di{z, fi,m) 
depends only on the mass and the factorization scale, but is insensitive to the high 
energy part of the process. In particular, the perturbative fragmentation function 
is universal, i.e. independent of the process. 

The PFF satisfy the DGLAP equations (2.29). The latter can be solved with 
NLL accuracy. To completely specify the solution however, one needs to specify at 
some scale //q an initial condition Dl^^{z, jiQ, m) that is also valid to NLL accuracy. 
Clearly, such a condition is also universal and can be obtained from a perturbative 
calculation [42] . To that end one first needs to observe that if the scale ji^ is chosen 
of the order of the mass m then no large logs will be present in the initial condition. 
Therefore, the initial condition can be simply calculated in perturbation theory: 




(2.58) 
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To obtain the above functions one needs to independently compute to order as 
both (7^ and a for some process and then plug the results into Eq.(2.57). Com- 
paring the terms order by order in as one gets in the MS scheme [42] : 



dl%(z,iio,m) = Cf 



1 + A. /^o 



ln^-21n(l-z)-l 



+ 



d^giz,^^o,m) = l(.^ + (l-.f)ln(^), 
d^X,gi',f,o,m) = 0. (2.59) 

Process independent derivations of these initial conditions also exist [43, 52]. 

The most convenient way to solve the DGLAP equations is to work with the 
Mellin moments of the fragmentation functions. Then the evolution equations 
take the following form: 
d 



,A(iV,/x,m) = ^^ij{N,as{pi))Dj{N,pi,m), (2.60) 



where we have introduced Di{N) through Eq.(2.39) and the anomalous dimensions 
are defined in Eq.(2.44). Also, one can conveniently factorize the dependence on 
the initial conditions: 

Di{N,i^,m) ^YlMN,fi,l^o) Df{N,i^o,m), (2.61) 
j 

where the functions Eij{N, /i, /iq) are universal solutions to the DGLAP equation 
(2.60) subject to the initial condition Eij{N, /iq, /iq) — Sij. Prom the factoriza- 
tion property (2.40) it follows that the 2;-space counterpart of this condition is 
Eij{z,no,no) = SijS{l - z). 

The DGLAP equations (2.60) couple all flavors that potentially participate 
in the process thus making the system hard to solve in general. An often made 
approximation is to neglect the mixing of the heavy flavor q with the light flavors 
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and the gluon. Such an approximation is justified in situations where the proba- 
bihty to create a heavy quark from a secondary gluon is smalL The approximation 
amounts to considering only the non-singlet contributions to the evolution kernels. 
We will consider only this case hereinafter. In this approximation it is very easy 
to solve Eq.(2.60) with NLL accuracy: 



p(i) 27r6i (0) 
On 



(2.62) 



Eq{N, /io) = exp I^P^t + [asii^o) - Q;s(/^)] 

where we have introduced the variable t as 

^^l^^aM_ (2.63) 

The coefficients 6o,i are defined in (2.15) and the functions Pjv'^^ are the Mellin 
transforms of the LO and NLO coefficients in (2.30). Their explicit form as well as 
iV-moments of the initial conditions (2.59) are given in the appendix. Collecting 
the resuhs in Eqns.(2.56), (2.57), (2.61), (2.62) and (A.27), we can write the 
differential cross-section for production of a heavy-flavored hadron as: 
1 da"" 



Q, m) = M^' C{N, Q, /x) Eg{N, /.o) /^f (iV, /xq, m)D^p{N) 

(7q uZ '- 

With A^~^ we denote the inverse Mellin transform: 



[Z). 

(2.64) 



1 pC+ioo 

f{z) = M-'[f{N)]{z) = — dNz-'^fiN), (2.65) 

Zm Jc-ioo 



which is typically performed numerically. We will give more details in the next 
Chapters where we discuss our applications. In Eq.(2.64) we have also introduced 
the iV-space coefficient function: 

UZ z 

(To Jo 

and the non-perturbative fragmentation function Dnp{N). The latter is extracted 
from fits to experimental data. We will present details for how that procedure 
works in the next Chapter where we consider the case of top decay. 



1 /"l drr'i 

C{N, g, //) = - / dz z^'-' — iz, g, //), (2.66) 
(Tfi "'0 az 
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Eq.(2.64) is the improved cross-section for production of a hadron H carrying 
longitudinal momentum fraction z, and with all quasi-coUinear logs ln(m^/(5^) 
resummed with NLL accuracy. As we will also see in the next Chapter in the 
example of top quark decay, the application of the resummation via the perturba- 
tive fragmentation function represents a serious improvement over the fixed order 
perturbative calculation. 

2.4.3 Soft-gluon Threshold Resummation 

Before we move to the discussion of our processes of interest, let us describe 
another, more subtle source of large logarithmic corrections that arise to all orders 
in perturbation theory and also need to be resummed. Those logs appear in 
processes containing massless and/or massive particles and are also of IR origin. 

Let us go back to Eq.(2.43). Since no masses were present there, wc concluded 
that such a renormalization group improved cross-section does not contain any 
large parameters and therefore can be unambiguously computed. We did not 
consider however its behavior as a function of the dimensionless parameter N, the 
Mellin conjugate of z. This parameter can be large since the limit of large z (i.e. 
z ^ 1) corresponds through the Mellin transform (2.39) to iV — >■ cxo. Therefore in 
the kinematical limit of large z one may expect terms ~ ln(A'") to appear in the 
conjugate space and to spoil the perturbative expansion in as- Indeed, a typical 
perturbative calculation of a hard cross-section shows terms like 

2n 

ag^c^ln-(7V). (2.67) 

m=0 

The terms with n < m < 2n are called leading logs (LL), the ones with m — n 
are next-to-leading (NLL) logs and so on ^. An example is the large behavior 
"^We will give more convenient characterization of these logs below; see Eq.(2.80). 
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of both space- and time-like anomalous dimensions shown in Eq.(2.36) where the 
latter are as singular as a single power of ln(A'"). Relations between the large z 
and large N behavior of various relevant functions are given in the appendix. 

The physical origin of the above mentioned (soft) logs is in the soft gluon 
radiation near the boundary of phase space. As we discussed in Section (2.4), the 
soft divergences cancel when all the contributions from real and virtual radiation 
are taken into account. However, this cancellation is incomplete near the boundary 
of the phase space 2; — > 1 for semi-inclusive cross-sections, and as a result large 
terms ~ ln'^(l — z)/{l — z) are present in the perturbative calculation. 

The all order resummation of classes of such soft logs is performed in terms 
of the Mellin moments because in Mellin space the factorization properties of 
the amphtudes are explicit. We will consider cross-sections of the type (2.64), 
where D can be a general fragmentation or distribution function. The form of 
Eq.(2.64) suggests the following approach for the resummation procedure: one 
can obtain separate resummed expressions for each of the factors appearing in the 
RHS of (2.64) that require such resummation. This way, one explicitly separates 
the process-dependent from the universal factors which makes the resummation 
procedure more universal. Of course, this is not the only possibility; one can 
solve directly for the whole cross-section. An example is the large-z resummation 
performed in [42]. 

In what follows, we will describe the method of Catani and Trentadue [53] 
originally apphed to Drell-Yan and DIS. We will apply their method in the next 
two chapters for the calculation of the resummed expressions for the coefficient 
functions for top-quark decay and inclusive CC DIS. The equivalence of their 
result with the Stcrman's [54] was shown in [55]. 

The coefficient function represents the normalized cross-section for production 
of a massless quark after the corresponding coUinear singularities are subtracted. 
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Up to non-singular terms in the limit ^ — >^ 1, it is approximated by the eikonal 
cross-section [53]: 

da^as^ [d<l>^^\n\S\0)\'', (2.68) 

n=0 

where / dcf)^"'^ is the (partially integrated) phase space for n-soft gluons, \n) is a 
state with n real soft gluons and S is the (nonlocal) eikonal operator: 



S^TPex.p 



-igs j dx^A^ 



(2.69) 



The contour C is the classical path of the radiating hard partons and A is the gluon 
field in the Heisenberg picture. T is a time-ordering operation for the operators A 
and P denotes path-ordering for the color matrices along to contour C. Prom the 
unitarity property ^ of the eikonal operators [56], it follows that the total eikonal 
cross-section equals exactly the leading order cross-section, i.e. the higher order 
contributions vanish (in particular that also implies cancellation of the soft-gluon 
divergences). For example at order as'- 



(0) 



+ cu(^)=0, (2.70) 



where uj^'''\ i = 0, 1 is the total probability for virtual {i = 0) and real (i = 1) soft 
gluon emission. 

One can evaluate the eikonal cross-section to some fixed order in as by work- 
ing in the eikonal approximation. In this approximation, the emission of a soft 
gluon with momentum q from a hard quark with momentum p is described by 

vertex 2p^ and with eikonal propagator l/2p.q. An important property of the 

eikonal approximation is that real gluon emission factorizes [40]. To order as the 

amplitude for a real soft-gluon emission in a process involving k hard particles is 
*For the applications we consider in this Thesis it is satisfied, since the unitarity of the 
eikonal operators is explicit in the case when every two points of the contour C are separated 
by a time-like interval, i.e. the contour is time-ordered; see [56] . 
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proportional to the /c-particle hard amphtude. In general however, there are color 
correlations involved. In the case k — 2,3 one can show that the correlations are 
absent since the product of color matrices can be expressed completely in terms 
of the Casimir (2.10) [40]. For the purposes of our apphcations we will only need 
the result for the case of two hard partons (quarks): 

where J'^ is the eikonal current: 

J^=pL_pL^ (2.72) 

Pi-q P2-q 

and pi_2 are the momenta of the two hard quarks that constitute the path C in 
Eq.(2.69). In case the hard partons are gluons the color factor in Eq.(2.71) should 
be replaced with Ca- 

Collecting the results above, one can easily obtain the eikonal coefficient func- 
tion to order as'- 

C^oft^z) = —^=(l + ^(0)) d{l -z)+ [ duj^'^S (l-z-^]+ 0{al), (2.73) 
gb dz ^ ^ J \ poj 

where qo is the energy of the soft gluon and po is the cm. energy of the hard 
parton. To obtain the all-order result valid up to NLL accuracy one has to take 
the Mellin moment of Eq.(2.73) and then to exponentiate the as term. The 
resulting integral typically takes the following form: 



+ J[a5((l-^)V)]}- (2.74) 

To put Eq.(2.74) in its present form, we have introduced the variable z — l — qo/po 
(as follows from the argument of the (5-function in (2.73)) and have worked in 
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spherical co-ordinates conveniently determined by the hard momenta. Two of 
the integrations over the phase-space of the radiated gluon are trivial. The one 
remaining integration is over the azimuthal angle which in turn can be related 
to the transverse momentum of the radiated gluon. The latter is known to set 
the scale of the running coupling in the soft region [57], [7] and is an important 
leading order effect. More details on these relations will be presented when we 
consider our applications. The origin of the —1 term in Eq.(2.74) is in the virtual 
emission term co^^^ after (2.70) is used. The number a is process dependent; 
typically a = 1,2. The kinematical constraints imply that the lower limit of the 
integration is zero. It is convenient to regularize the low-g^ collinear divergence 
by introducing a cut-off. The (singular) dependence on the cut-off is later absorbed 
into the partonic distribution/fragmentation function as was discussed in Section 
(2.4.1). As a result, in the MS factorization scheme the lower limit is set by the 
factorization scale fip [43], [58]. 

The function J is process dependent and has a perturbative expansion in as- 
To NLL accuracy, only the 0{as) term participates. That function receives two 
contributions: the first one is from the coUinearly non-enhanced term of the square 
of the eikonal current, after the integration has been explicitly performed. Up 
to NLL accuracy one can use the following result: 



To show Eq.(2.75), one needs to use the running of the strong couphng (see 
Eq.(2.17)) and then to carefully identify and collect the NLL contributions. Both 
in (2.75) and (2.74), Q stands for some process-dependent hard scale. 

To specify the second contribution to the function J, we need to first explain 
the physical origin of the various enhanced contributions in the soft limit. From 
kinematical considerations it is clear that in the limit z — > 1, typically, the initial 




a 



(2.75) 
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partons can emit only soft gluons. The hard partons in the final state however 
may or may not be subject to that constraint. In the case of inclusive DIS for 
example, the hard parton in the final state is completely free to fragment with 
the only restriction being on the virtuality of the final jet. If however the final 
parton is observed (and considered as massless as is appropriate for a coefficient 
function) it can only radiate soft gluons. Since the eikonal approximation handles 
only the contribution from soft gluon radiation (that itself can be coUinear), one 
needs to specify in addition to the eikonal cross-section the contribution from the 
fragmentation region. Such a contribution was evaluated in [53] after the soft- 
gluon radiation was completely decoupled from the final massless quark (with 
appropriate choice of the vector specifying the physical gauge). The jet function 
that incorporates the effect of the coUinear radiation in the fragmenting jet was 
obtained from a modified evolution equation. We will return to that point later 
when we consider our applications. 

Finally, the results to NLL order can be completed by introducing an additional 
(9(q;|) contribution to the RHS of Eq.(2.74). One needs to make the replacement: 



TT q"^ q^ 



(2.76) 



where the function A{pLs) can be expanded as follows: 

A{as)^f^['^\ A^-\ (2.77) 



n=l 



TT 



The first two coefficients are needed at NLL level and are given by [53, 59]: 

A(^) = Cf, (2.78) 

(2.79) 



2 



I 18 6 / 9 ^ 



A numerical estimate for A^^^ is known too [60] (see also [61]). The function A 
represents the coefficient of the singular term of the anomalous dimension (2.36). 
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After the replacement (2.76), Eq.(2.74) contains all NLL contributions. Its 
explicit evaluation would, however, contain subleading terms as well. One can 
evaluate explicitly the integrals and extract only the LL and NLL terms. After 
that is performed, the eikonal cross-section takes the form: 

C'"f\N, Q, as{fi'),fi, fiF) = exp [in Ng^'\X) + y(2)(A, Q, fip)] (2.80) 

where we have specified all arguments, restored the separate dependence on the 
renormalization and factorization scales, and introduced: 

X^boas{iJ,'^)lnN. (2.81) 

The functions g^^^ and y^^) contain all LL {]n{N){as In(Ar))^) and NLL {{as In(Ar))'^) 
contributions respectively. The results at NNLL order and beyond can be also 
written in terms of functions Q;|(A^)g'('^+^)(A), k — 1, . . . (see e.g. [62]). To eval- 
uate the integrals in (2.74) with NLL accuracy, one makes use of the running 
coupling Eq.(2.17) and the following identity, valid to NLL [53]: 

.AT-i _ 1 ^ _e (^1 _ £^ _ . (2.82) 

We would like to emphasize the fact that the use of the replacement (2.82) makes 
it possible to evaluate the integral in (2.74). It is clear that the integrals in (2.74) 
are divergent for any value of N simply because of the presence of the Landau pole 
in the running couphng as; see Eq.(2.16). The role of the replacement (2.82) is 
that it extracts finite values from the integrals in (2.74) up to NLL accuracy, while 
neglecting subleading terms that are formally divergent. Consistency is however, 
preserved because of the neglecting of similar subleading terms in the derivation 
of (2.74). A detailed description can be found in [63] and [64]. 

So far we discussed only the evaluation with NLL accuracy of the coefficient 
function. However, the 0{as) expression for the initial condition of the quark 
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initiated perturbative fragmentation function (2.59) also shows enhancement in 
the hmit 2; — > 1. The resummation of the soft-gluon contribution to the (universal) 
initial condition for the PFF was performed in [43] with NLL accuracy. We shall 
only present the final result here: 



+ H 



1-z 
as ((1 - zfm^) 





(l-^)2m2 



(2.83) 



Here m is the mass of the heavy quark and the function H also has a perturbative 
expansion in. ag'- 

H{as) = Y^('^X H^'^y (2.84) 



n=l 



TT 



At NLL level, we are just interested in the first term of the above expansion: 



(2.85) 



The explicit evaluation of the initial condition (2.83) to order NLL gives: 



LT^exp ln7V52(Ao) + <7ln«(Ao) 



.(2)/ 



(2.86) 



where we have suppressed the arguments for brevity, and Aq is given in (2.81) with 
the replacement //^ ^ /Xq- The functions g^f^ are: 



2nboXo 
AW 



[2Ao + (1 - 2Ao) ln(l - 2Ao)] , 



+ 
+ 



27r6o 
1 

27r6o 
27r6o 



In ^ + 2jE 



ln(l - 2Ao) 



4Ao + 2 ln(l - 2Ao) + ln2(l - 2Ao) 



[2Ao + ln(l - 2Ao) 



ln(l - 2Ao) 



+ In 



(2.87) 



(2.88) 
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So far we discussed how one can evaluate the approximate expressions for 
the coefficient function and the initial condition of the PFF in the limit z — > 
1. However, to obtain a result equally applicable for all z, one needs to give a 
matching prescription for the fixed order result and its approximation in the soft 
limit. In this Thesis we will use the following matching procedure: one adds the 
fixed order result and the result valid in the soft limit together, and then subtracts 
what they have in common. Note that for consistency, the decomposition of the 
result valid in the soft limit in powers of as should coincide with the soft limit of 
the fixed order result. We will give details for that procedure in the next Chapters. 



Chapter 3 

6-quark Fragmentation in t-quark 
Decay 

In this Chapter wc present our original results for ^-fragmentation in top decay. 
In their derivation we apply several techniques within perturbative QCD such as 
the fixed order calculation and resummation of large quasi-coUinear and soft logs. 
Those techniques were discussed in the previous Chapter. We also extract the 
non-perturbative fragmentation function from fits to e'^e~ data and thus make a 
prediction for the spectrum of 6-fiavored hadrons in top decay. This Chapter is 
based on the papers [51] and [65]. 



3.1 Motivation: Top Physics 

The existence of the top quark is a prediction of the Standard Model. It was 
discovered at the Fermilab Tevatron [66] and its mass was measured [67] to be 
174.3 ±5.1 GeV. This large mass applies, among other things, a large coupling 
to the Higgs boson (the top- Yukawa couphng quh ~ 1)- Unhke the other five 
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quarks, top's CKM matrix elements are dominated by a single element ~ Vtb - 
which is very close to one. As a result top decays almost exclusively to a real W 
boson and a 6-quark. Its electric charge (as predicted by the SM) is +(2/3)e but 
is not yet measured directly [68] . The present status of these and other top-quark 
properties can be found in the extensive recent review [69]. 

At present, and for the foreseeable future, the study of processes involving 
the top quark is of great interest. On the experimental side, because of its very 
large mass, the top is still a challenge to present experimental facihties. The CDF 
and D0 collaborations in Run I at the Fermilab Tevatron produced a small but 
sufficient number of top events to directly prove its existence in 1995 [66] and 
to pinpoint its mass. Since then the Tevatron has been upgraded and Run II is 
presently collecting data. The expectations are that during Run II more detailed 
measurements of the top quark properties will be possible. In the near future 
(around 2007-8) the LHC proton-proton coUider is expected to start operation. 
Due to its much larger energy of operation it is expected to produce vast numbers 
of top quarks with a corresponding improvement in top studies [70]. There is 
another high-energy machine - the electron-positron Linear Collider [71] - the 
prospects of building which are being studied at present. Such a machine will 
not produce as many top events as LHC, but will have the potential for precise 
measurements of the top quark parameters. 

On the theoretical side, the large mass of the top quark makes it decay very 
fast to a 6-quark and a boson. That imphes a very large top decay width, 
the NLO value being Tfop ~ 1.4 GeV. The corresponding hfe-time is about an 
order of magnitude shorter than the typical hadronization times set by the scale 
Aqcd- An immediate consequence is that the top quark cannot form t-flavorcd 
hadrons; it simply decays before it is able to hadronize [72]. For that reason the 
top behaves much hke a free particle. That property of the top as well as the very 
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small ratio of its width to its mass suggest that it is reasonable for most purposes 
to assume that top production and top decay factorize to a good approximation. 

It is not only the study of the SM that draws attention to the top quark. As 
we mentioned, the top has very large Yukawa coupling, i.e. its coupling to the 
Higgs sector is largest among all known particles. For that reason it is reasonable 
to expect that new physics that may be related to the mechanism of electroweak 
symmetry breaking may first show up in the top sector. For a more detailed 
account of such models see [69] . 

Clearly, a precise knowledge of the properties of the top quark is necessary 
not only in order to verify experimentally the SM predictions but also to be 
able to discriminate non-SM effects should they be discovered in forthcoming 
experiments. However, because the top has such a small life-time, the only way 
one can get any information about the top itself is through its decay products. 
Clearly, this requires - among other observables - a reliable description of the 
fragmentation of the 6-quark in top-quark decay. 

As shown in [73], 6- fragmentation is indeed one of the sources of uncertainty 
in the measurement of the top mass at the Tevatron. At the LHC, recent studies 
[74] have suggested that final states with leptons and a J/^', with the J/'^ coming 
from the decay of a 6-flavored hadron and the isolated lepton from the W decay, 
will be a promising channel to reconstruct the top mass. In [74], the expected 
experimental error, including statistics and systematics, has been estimated to 
be Arrit ~ 1 GeV and the b fragmentation is the largest source of uncertainty, 
accounting for about 0.6 GeV. 

In this Chapter we will study the 6-spectrum in top decay in the following 
framework: neglecting the interference effects between top production and decay. 
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one has 



1 da 



1 dV 

r dxb 



(3.1) 



a dx, 



b 



In the above equation, {l/a)da/dxh is the normahzed differential cross section 
for the process of production of a 6 quark with energy fraction Xh (defined in 
Eq.(3.5)) via top quarks for any production mechanism. Similarly, {l/T)dT / dx^ is 
the normalized differential distribution of 6-quark in the decay of a real top quark. 
Our purpose in this Chapter is to give a precise prediction for the distribution on 
the RHS of Eq.(3.1); through the relation (3.1) our results will then be applicable 
to top quarks produced in pp (Tevatron), pp (LHC) or e^e~ (Linear Collider) 
collisions. 

We start by first calculating in Section (3.2) the differential top width with 
respect to the 6-energy fraction in NLO QCD taking full account of the mass 
of the h quark. Then in Section (3.3) wc resum with NLL accuracy the large 
quasi-collinear logs ln(m^/m|). Applying soft-gluon resummation techniques in 
Section (3.4) we improve the prediction for the 6-energy spectrum at large — > 1. 
Then we discuss the improvement to the fixed order calculation due to the various 
resummations. Finally, in Section (3.5) we extract the non-perturbative fragmen- 
tation function using e^e~ data to make predictions for the energy spectrum of 
6-flavored hadrons in top decay. 



In this section we present our results for the NLO QCD calculation of the differ- 
ential width: 



3.2 t^hW \Yi NLO QCD 



1 dVf, 
Lo dxb ' 



(3.2) 
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for the Standard Model top-quark decay: 

t{q, Tilt) b{pb, mb)W{pw, mw) {g{Pg: 0)) , (3.3) 

of an on-shell top quark at next-to- leading order in as, with respect to the en- 
ergy fraction of the produced 6-quark (defined below). We assume that the de- 
cay process (3.3) is the only decay mode of the top quark. That assumption 
is consistent with recent measurements of the CDF Collaboration of the ratio 
R — B(t — > Wb)/B(t — > Wq), where q is a d, s or b quark, and the subsequent ex- 
traction of the Cabibbo-Kobayashi-Maskawa matrix element Vtb [75]. We present 
the results from two independent calculations: when the exact dependence on 
the bottom mass is kept, thus regulating the coUinear divergences, and when the 
bottom mass is set to zero from the outset. In the latter case additional coUinear 
divergences occur. In all our calculations we use dimensional regularization to 
regulate the UV and all IR divergences, that is, we perform the calculations in 
D = 4 — 2e space-time dimensions (see appendix) . We use the on-shell renormal- 
ization scheme to remove the UV divergences. Let us note that in principle one 
does not need to calculate Eq.(3.2) both in the nib — and rrib cases; one can 
infer the massless result for the normalized distribution (3.2) (with the coUinear 
singularity subtracted in the MS scheme) from the massive one, using the initial 
conditions for the PFF, Eq.(2.59). In fact with our independent calculation of 
those two cases we were able to reproduce those initial conditions in the case of 
top-quark decay. 
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3.2.1 Calculation with m^, 7^ 



In this section we will present in more detail the calculation in the case when the 
full dependence on the 6- mass is kept. Let us first introduce our notation: 

b = ^ 



m 



w — 



s — 



w 



Q 

Go 



1 1 
mf 

lil + b-w), 

Vb 

s 

1 



1 + 6 - 2w + 



xl — P"^ — arcth 




where: 



Xe 



Pb-Pt 
mfs 



P<Xb<l. 



(3.4) 



(3.5) 



Xe ,max 

It is obvious that Xb is an invariant, which in the top rest frame coincides with 
the energy of the 6-quark, normalized with respect to its maximum value. The 
LO width To is: 



amt\Vtb\ 



4gG'o, 



(3.6) 



16 sin^ 9w 

where a is the electromagnetic coupling constant and 6w is the Weinberg angle. 

Next we proceed with the details of the calculation. At order 0{as) one has 
to calculate the contributions from the diagrams in Figures (3-1), (3-2) and (3-3). 
To obtain the contribution from the real emission diagrams - Fig. (3-3) - one has 

to add the two matrix elements and square them. After averaging/summing over 
spins and colors (all the needed details can be found in the appendix) one gets 
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Figure 3-1: Tree level (Born) diagram. 
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Figure 3-2: Contributing virtual emission diagrams. 



< 



W 




g 




Figure 3-3: Contributing real emission diagrams. 



the following expression for the square of the matrix element: 



\Mreal\^ 



2 2 iTA |2 

3 5 9s\Vtb\ 



III 



_{Pt-Pg){Pb-Pg) 



{2sGo - e4s^ 
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+ 



1 



l + b, ,2 1 + 6' 



w 



-{pt.pgy + [2 + 



w 



{Pt-Pg){Pb-Pg) 



1 



+ 



(Pb-PgY 



^ ^ ^{Pb-Pgf + ( 2 + ) {PfPg)iPb-Pg) 



w 



w 



2{pb.pg)m'^Go + m^i-bGo + e2bs) 



(3.7) 



in terms of the scalar products of the four vectors. We have kept those terms 
hnear in e that will give finite contributions after the singular integration over the 
phase space is performed. The next step is to integrate the matrix element (3.7) 
over the three-body phase space. More details about that procedure are relegated 
to the appendix. As a result, we obtain the following contribution from the real 
emission diagrams to the differential width (3.2) at order as'- 



1 dV 



real 



To dxb 



ttQ 

+ s(Li2| 



<^{xb^l) (--27E + 2 



In 



(Q 



2Q 



1-s + Q^ 



— Lio 



2Q 



- In 
+ In 



. Vb . 



w 



2 Hxb) ^ 



s-6 + 2^j +ln 



s-b+Qj 

1-s + Q' 



2s{l-P)j 



s 



+ 



1 + 



(5(1 - Xb) 
1 + 



(3. 



+ 



2V^6 - 



■ Xb). 

2 — 



Go 

1-Xb 



2w 



1- Xb) -1 

1 + 1 



(1 -Xb)-1 



Go \l-2sxb + bj ' Gn\~ ' 2w 

As expected, the real emission contribution to the differential distribution 
contains an explicit soft divergence ~ 1/e. The latter will be cancelled after we 
add the contribution from the virtual diagrams. Let us explain their evaluation. 
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The diagrams with virtual gluon emission that contribute to order as are 
shown in Fig. (3-2). They represent the same physical process as the leading order 
(Born) diagram in Fig. (3-1). Because the virtual diagrams contain loop integrals, 
they are UV divergent. We regularize these divergences by working in D space- 
time dimensions. The cancellation of the UV divergences requires an appropriate 
introduction of counterterms. The counterterms are determined from the on-shell 
condition: 

E(^=m) = 0, ' ^ 0, (3.9) 

where S is the quark self-energy for quark q with mass m. From the second con- 
dition in (3.9), one can derive the quark wave function renormalization constant 
(it coincides with the one in [76]): 



This then determines the vertex renormalization constant: 

S7r\euv eiR 2 {ATr/iy 2 \47riiy J ^ ' 

The renormalized vertex thus becomes: 

(Z(^) - l) 7^. (3.12) 

The function contains the higher order contributions to the vertex amplitude, 
and is defined through the relation: M — g / \/9>VthuK^{l — ^^)ue'^. 

Let us make a remark on the evaluation of the second term in Eq.(3.9). The 
function ^(p) is UV divergent and its finite part is a regular function of the 
momentum p for < m^. The point p = m is a (finite) logarithmic branching 
point, i.e. there E contains terms behaving as: 

lini E ~ y In (y) -|- regular terms. 
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For that reason, although E is finite for p = m, its derivative with respect to p 
diverges at that point. This is an IR divergence which we have regulated and 
distinctly parameterized in Eq.(3.10). 

The evaluation of the vertex correction and of the self energy diagrams can 
be done by applying the Passarino-Veltman method [77] for reduction of tensor 
integrals. As a result, after some extensive algebra, one can reduce the calculation 
to the evaluation of the three scalar integrals Aq, Bq and Cq. The latter are 
defined as the integration in D space-time dimensions of a product of one, two 
and three scalar propagators with, generally, different masses. For their evaluation 
we use standard techniques and Feynman parametrization. In the present case of 
one-to-two particle decay, the phase space is very simple (see the appendix). To 
obtain the differential cross-section one does not need to perform any nontrivial 
integrations. The resulting expression is: 



Fo dxb ttQ [ \e \47r^/x^ 

+ s 



1-s + gy \s-b + Q^ 

l-s + Q\ , fs + Q' 



2 m (. + g) m j - m m ((. + Q) v.; 

- 2s' , /s + Q\ fs 



Q 



AwC 



-(1 + 46 + 6w - 56^ -5w + Aw'^) In (Vb) \s{l - Xb) (3.13) 
J 



To get the needed distribution we need to only add the contributions (3.13) 
and (3.8). After also adding the LO contribution, the result reads: 



Fo dxb ttQ 



Lio I - — — — -I - Lio '^^ 



1-s + Ql ^\s-b + Q^ 
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+ lln(6)lnfi-A±|^ 



+ (1-6) In 



w 



+ Q 



6- 



[w 



-bj\n\^ 
b){s-b) 



Vb 



wGo 



- 1 



ln{b) 



(3.14) 



- 2 In 
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2 $(xb) 



(1 - Xb)- 



2t 



+ 



5(1 -Xb) 

1 + 6' 



1 + 



Go 
1 - 



2W 



l-Xb)-l 
1+6 



2«; 



;i-a;6)-l 



Go Vl - 2sxb + 6. 

Separately, the contributions from the virtual and from the real emission are IR 
divergent but those divergences cancel in the sum. In the derivation of the above 
result, we have extensively made use of various relations between the dilogarithms 
(Spence functions); see the appendix for details. Finally, one can verify that the 
integral of the differential distribution (3.14) coincides with the total top-decay 
width to order as evaluated in [76]. As expected from the KLN theorem (see 
Section (1.4)), all IR sensitive contributions disappear. Let us also note that the 
top width is known at present to order a| [78]. 

Let us also present here the result (3.14) in the limit 0, where all power 

corrections of the type mf,/mt can be neglected. The result reads: 



= d(l - Xb) + — — Ai{xb), 



To dxb 



with 

Ai{xb) 



= C, 



27r 



mi 



(3.15) 
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(3.16) 
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The origin of the plus-prescription appearing in the above results as well as some 
its properties are discussed in the appendix. 



3.2.2 Calculation with mb=0 



In a similar fashion to the above derivation one can also obtain the massless 
results with the bottom mass set to zero from the outset. In this calculation we 
proceeded in the same way as in the ruf, ^ case. The integrals to be evaluated 
simplify significantly in this case. We use the same renormalization scheme as 
in the case nib 0- The essential difference is that in the massless case, the 
contributions to the differential width from real and virtual gluon emission, after 
the UV renormalization is performed, contain IR divergent terms up to 
An extra power of 1/e appears because the 6-mass no longer prevents coUinear 
divergences. For that reason, in general, one need to keep terms to order O(e^). 
Similarly to the case ^ 0, in the sum of the real and virtual contributions 
to the differential width the leading IR divergences ~ cancel. The resulting 
expression contains a singularity ~ 1/e which is in one-to-one correspondence with 
the quasi-coUinear log ln{mb/mt) present in the massive result (3.16). Finally, the 
massless result reads: 

' S{1 - -.) + ¥^ ic, + ^^(1 - xb) 



Fo dxb 



27r \ 



X 



-- + 7i5-ln47r) +Ai{xb)}, 



(3.17) 



with Ai (xb) defined as: 
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4z 



+ 2(1 + ^2) 



1 - 



w{l - w){l - zf 



(1 + 2w)(l - {l-w)z) 



+ S{l-z) 
2w 



4Li2(l -w) + 2 ln(l -w)\nw-^ + llL^ in(i _ 



1 3w 
m w + ^ z 9 



(3.18) 



1-w l + 2w 

In order to get the correct finite term in the normahzed differential width, the Born 
width To (see Eq.(3.6)) will have to be evaluated in dimensional regularization to 
order 0{e). We find: 

amt\Vtb\^ {l-w)\l + 2w) 



To 



IGsin^ ew 



w 



X 



-7e + In 47r - 2 ln(l - w) + 2 



{1 + e 

1 + w 



l + 2w 



(3.19) 



One can easily verify that the difference Ai{xh) — Ai{xh) exactly reproduces 
the (heavy quark initiated) initial condition for the perturbative fragmentation 
function, Eq.(2.59). We also checked that the integral of Eq. (3.17) agrees with the 
result of [79], where the 0{as) corrections to the top width have been evaluated 
in the approximation of a massless b quark. Of course, it also coincides with 
the massless limit of the total width obtained from our calculation with non- 
zero bottom mass. It is worth noting that in the total width the IR divergences 
disappear together with the dependence on the unphysical factorization scale /ip- 
Technically that is obvious since the terms that depend on hf are multiplied in 
(3.17) by the LO Altarelli-Parisi splitting function (2.34) whose integral vanishes 
(see Eq.(2.31)). 
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3.3 NLL Resummation of logs ln(m^/m^) 

In principle the perturbative result Eq.(3.14) gives a finite prediction for the b- 
spectrum in top decay. However, that quantity is not IR safe (in the sense of 
Section (2.4)) since it diverges in the limit — > 0. As follows from our general 
discussions in the previous Chapter, to obtain a reliable perturbative prediction 
one needs to resum the quasi-collinear logs \n{ml/mf) to all orders ^ in the strong 
coupling. We will show later in this section that the numerical difference between 
the two expressions for the 6-spectrum - Eqns.(3.14) and (3.15) - is totally neg- 
ligible. Since the two expressions differ by power corrections ~ (mft/m^)" which 
are small since mb/rrit — C(10~^), one can conclude that the power corrections 
can indeed be safely neglected. Let us recall that this is precisely the condition 
for validity of the method of the perturbative fragmentation function described in 
detail in Section (2.4.2). Therefore, in the following we will present the applica- 
tion of this method for resummation of the above mentioned quasi-collinear logs 
in top-quark decay with NLL accuracy. 

Following the general discussion in Section (2.4.2), the differential width for 
the production of a massive b quark in top decay can be expressed via the following 
convolution: 

1 dr 1 f^dzdr M^f^b \ , X 

— — (X6,mt,mw^,m6) = — / — — {z,mt,mw,IJ',IJ'F)D[—,iiF,mb), (3.20) 
To dxb To Jxb z dz \ z J 

where we have restricted the above result to the flavor non-singlet case, i.e. to 
the perturbative fragmentation of a massless b into a massive b. In doing so, we 
have ignored the strongly suppressed mixing of the 6-quark with the gluon and 
the other flavors through DGLAP evolution. In Eq. (3.20) fi and /ip are the 
^ There is another mass parameter in our problem - the W mass mw- However, we will treat 
it as not small with respect to mt, and therefore will not resum terms like ln(TO^/m|). 
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renormalization and factorization scales respectively, the former associated with 
the renormalization of the strong coupling constant. In principle, one can use two 
different values for the factorization and renormalization scales; however, a choice 
often made consists of setting fi = /ip and we shall adopt this convention for most 
of the results which we shall show. 

The NLO, MS-subtracted coefficient function is: 

- MS 

'^"^^ S{l-z) + ^Mz), (3.21) 



To dz ' ' 27r 

with ii given in Eq.(3.18). The consistent NLO plus NLL (as defined in Eq.(2.51)) 
evaluation of the 6-production rate requires that «5 is the NLO strong coupling 
constant (2.16) with Uf — 5 for the number of active fiavors. The initial condition 
for the perturbative fragmentation function must be evaluated to order as as well. 
The initial condition is presented in (2.59) with m being the bottom mass and, 
in order to avoid large logs there, we need to take the initial scale /iof of the order 
of nih- 

As follows from our discussion of Eq.(2.64), we will work in analytical form in 
terms of the Mellin moments N; the Mellin transformation is defined in (2.39). 
The Mellin transform of Eq. (3.20) is: 

TNimt, niw, rub) = VN{mt, mw, A*, IJ'F)Db,N{lJ'F, mb), (3.22) 

with 

Vj^{mt,mw,rnb) ^ ^ I dxb x^~'^^{xb,mt,mw,rnb). (3.23) 
i -^0 aXb 

The fragmentation function Db^^ is a product of the solution of the DGLAP 
equations Eq.(2.62) with the initial condition Eq.(A.27). The exphcit A^-space 
expression for the coefficient function f jv can also be found in the appendix. 

For our numerical study, we shall assume rrit — 175 GeV, mw — 80 GeV, 
mfe = 5 GeV and A = 200 MeV. 
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The 6-quark energy distribution in Xfe-space will finally be obtained by inverting 
the A^-space result (3.22) numerically, by a contour integral in the complex N- 
plane. We will use the so-called minimal prescription introduced in [63], where 
the integration contour is deformed with respect to the one in Eq.(2.65). We will 
discuss that point in more detail in the next section. 

We already have all the ingredients needed to obtain the expression for the 
6-quark spectrum in top-quark decay with coUinear logs ~ In (mt/mt) resummed 
with NLL accuracy. We normahze our plots to the total NLO width F, obtained 
by neglecting powers ~ {rrih/mty. Its expression can be found in [79]. In fact, 
that normalization is not affected by the factorization on the right-hand side of 
Eq. (3.20), the DGLAP evolution of the relevant factors, or the resummation of 
soft logarithms in the initial condition of the perturbative fragmentation function 
and the coefficient function. 

We first present our results for the Xb spectrum in Fig. 3-4. We plot the 
Xb distribution according to the perturbative fragmentation approach. For the 
sake of comparison, we also show the exact 0{as) result for a massive b quark, 
Eq.(3.14). We set /i — fxp — rrit and /iq — /iof — iTib- As an illustration, the 
effect of NLL soft-gluon resummation only in the initial condition of the pertur- 
bative fragmentation function is also presented; it will be systematically discussed 
in the next section. Clearly, the use of perturbative fragmentation functions has 
a strong impact on the Xb distribution. The fixed-order result lies well below 
the perturbative fragmentation results for about 0.1 <Xb< 0.9 and diverges once 
X(, — > 1, due to a behavior ~ 1/(1 — Xb)+- Moreover, as we previously noted, the 
full inclusion of powers of nib/mt has a negligible effect on the Xb spectrum; the 
dot-dashed and dotted lines in Fig. 3-4 are in fact almost indistinguishable. As 
for the perturbative fragmentation results, the distribution shows a very sharp 
peak, though finite, once Xb approaches unity. This behavior will be smoothed 
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Figure 3-4: 6-quark energy distribution in top decay according to the perturbative fragmenta- 
tion approach, with (sohd hnc) and without (dashes) NLL soft-ghion rcsumrnation in the initial 
condition of the perturbative fragmentation function, and according to the exact NLO calcula- 
tion, with (dot-dashes) and without (dots) inclusion of powers oi rrib/mt- In the inset figure, we 
show the same curves on a logarithmic scale. 

out after we resum the soft NLL logarithms appearing in the initial condition of 
the perturbative fragmentation function and the coefficient function; an indication 
is the behavior of the distribution after the soft-resummation (Eq.(2.86)) in the 
initial condition (2.59) is performed. Both perturbative fragmentation distribu- 
tions become negative for Xb ^ and xi, ^ 1, which is a known result, already 
found for heavy-quark production in e'^e~ annihilation [43], [46]. For Xb — > 0, the 
coefficient function (3.21) contains large logarithms ~ cks{iJ') Inxb which have not 
been resummed yet. Likewise, in the soft limit Xb — > 1, Eq. (3.18) contains contri- 
butions ~ as{ij)/{i — Xb)+ and ~ as{ij)[ln{l — Xb)/{1 — Xb)]+ whose resummation 
is discussed in the next section. As stated in [43], once Xb gets closer to unity, 
non-perturbative contributions also become important and should be taken into 
account. The region of reliability of the perturbative calculation at large Xb may 
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Figure 3-5: (a): a;6 spectrum for /xf = and /Uqf = mb/2 (dots), //of = nib (solid) and 
/iOF = 2m6 (dashes); (b): /xof = and /xf = mt/2 (dots), = mj (solid) and fip = 2mt 
(dashes). The renormalization scales are kept a.t n = rrit and no = nib. Though not visible, all 
distributions show a finite, sharp peak once Xb is close to 1. 

be related to the Landau pole in the expression for the strong coupling constant 
(2.16), and is estimated to be < 1 — A/mi,. Fig. 3-5 shows the dependence of 
the perturbative fragmentation Xb distributions on the factorization scales /iof (a) 
and fip (b). It is obvious that the results exhibit non-trivial dependence on the 
choice of those non-physical scales. Still, this dependence is concentrated close 
to the region of the Sudakov peak and is relatively mild far from it. Since the 
presence of scale dependence indicates the importance of higher order terms that 
are unaccounted for, we expect that after we resum large soft-logs the result will 
be less sensitive to the choice of those scales. That will be discussed in the next 
section. As a check of our results we have verified that for hf approaching /^qf, 
the distribution gets closer to the fixed-order, unevolved one shown in Fig. 3-4. 
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3.4 NLL Threshold Resummation 

The MS coefficient function (3.21) and tlie initial condition of the perturbative 
fragmentation function (2.59) contain terms behaving like 1/(1 — Xft)+ or [ln(l — 
Xb)/{1 — Xb)]+, which become arbitrarily large when approaches one. This is 
equivalent to contributions proportional to In N and In^ N in moment space, as 
can be seen by writing the MS coefficient function in the large- limit ^: 



ml 



+ 2 - 41n(l - w) - 21n ^ 



In AT 



+ X(mt,mi^,/ii.) + oQ) I (3.24) 

In Eq.(3.24) we have introduced the function KlmtjUiw, /J'f), which contains 
terms which are constant with respect to N; more details for the extraction of the 
terms which ~ const when N ^ oo can be found in the appendix. The function 
K reads: 

(3 \ 171^ 

- - 27E In ^ + 27I + 27E [1 - 21n(l - w)] 

1 — w 'Zw 

+ 21nu'ln(l - w) - 2 In(l-'u;) Inw 

^ ' l + 2w ^ ' 1-w 

+ 4Li2(l - w) - 6 . (3.25) 

3 

The Xh ^ 1 {N ^ 00) limit corresponds to soft-gluon radiation in top decay. 
These soft logarithms need to be resummed to all orders in as [53] , [54] to improve 
our prediction. 

We can perform the soft-gluon (threshold) resummation for the coefficient 
function following the general procedure described in Section (1.4.3). We treat 
the 6-quark as exactly massless and therefore the result is coUinearly divergent. 

^Following [62], we note that, by defining n = Nexp{'yE), we could rewrite this expression 
in terms of ln(n) rather than In N, with no 'Je terms explicitly appearing. 
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That divergence is subtracted in the MS scheme. The eikonal current reads: 

rrit „ {Pt ■ Pb) 



\J{.PuPb,Pg)Y 



-2- 



(3.26) 



{PfPgf (Pt ■ Pg){Pb ■ Pg) 

We express the 0{as) width in the soft approximation as an integral over the 
variables = (pb + PgYxg and z = 1 — Xg, with < < — — z)"^ and 

< ^ < 1. The limits z ^ 1 and ^ correspond to soft and coUinear emission 
respectively. In the soft approximation, z x^] Xg is the gluon energy fraction 
that is defined similarly to X}, (3.5). In the m^ — Q case which we consider here: 

1 2pfe • Pt I 2pg ■ Pt 

Xb = 2~' ^9 = ~^ 2~' < Xb,g < 1. (3.27) 

We point out that our definition of the integration variable is analogous to 
the quantity {1 — z)k'^ to which the authors of Ref. [53] set the scale for ag for 
soft-gluon resummation in Drell-Yan and Deep- Inelastic-Scattering processes. For 
small-angle radiation, ~ q^, the gluon transverse momentum with respect to 
the 6-quark line. The variable z is analogous to 2; = 1 — Eg/ Eg of Ref. [53] . 
Performing the operations discussed in Section (1.4.3) we obtain: 



tN{mt,mw:lJiF) = — / dz— [ 

TT Jo 1 — Z [Jfj.^ 

I /■m2(l-U))2(l-2'l2 



^m?(l-«;)2(l-z)2 dq 



2 



2 



/ dq as 

Jo 



(3.28) 



mf(l - wy{l - zf 

We note that it corresponds to Eq.(2.74) with a = 2 and = m^(l — w)"^. The 
origin of the second term in (3.28) is in the coUinearly non-enhanced term in the 
squared eikonal current (3.26). 

In order to include all NLL contributions we only need to perform the re- 
placement (2.76) and use the property (2.75). No additional contribution from 
coUinear radiation from the final state is presented here because of the character 
of the cross-section that our coefficient function represents; the final 6-quark is 
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"observed" and has its virtuality fixed by the on-shell condition pi = 0. Expo- 
nentiating the result we obtain the following form-factor that describes the soft 
limit of the corresponding production rate: 



,7V-1 



In A 



dz- 



1 r .m2(l-«,)2(l-z)2 ^^2 



(3.29) 



(3.30) 



^0 

+ S[as{m^ti^-'^fi^- 4' 
The function S{as) can be expanded according to: 

At NLL level, we are just interested in the first term of the above expansion: 

^(1) = -Cp. (3.31) 

The integral in Eq.(3.29) can be performed, up to NLL accuracy, according to the 
methods described in Section (1.4.3). This leads to the following result: 

AN{mt,mw,cxs{fJ^^), fJ^, fJ^p) = exp lnNg^^\X) + g^^\\, fj,, ^p) 



with A defined in (2.81) and the functions g'^^^ and g^"^^ given by: 



+ 



27r6oA 
AW 

27r6o 
1 



[2A + (l-2A)ln(l-2A)] , 



(3.32) 



(3.33) 



In -^^-^ — - 27E 



ln(l - 2A) 



4A + 2 ln(l - 2A) + ln^(l - 2A) 

^(2) 



+ 



27r6n 



[2A + ln(l - 2A)] 
ln(l-2A). 



(3.34) 



In Eq.(3.32) the term \a.N g^^^{X) accounts for the resummation of leading log- 
arithms In**"*"^ in the Sudakov exponent, while the function g^'^\X, fj,, fxp) 



resums NLL terms In" N. 
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Furthermore, we follow Ref. [43] and in our final Sudakov-resummed coefficient 
function we also include the constant terms of Eq.(3.25): 

asifJ,^) Cf 



1 + 



27r 



-K{mt,mw:I^F) 



X exp[ln7Vg«(A)+^(2)(A,//,//p)]. (3.35) 



One can check that the 0{as) expansion of Eq.(3.35) yields Eq.(3.24) as it must 
for consistency. Next we match the resummed coefficient function to the exact 
first-order result, so that also 1/N suppressed terms, which are important in the 
region Xb < I, are taken into account. We adopt the same matching prescription 
as in [43]: we add the resummed result to the exact coefficient function and, 
in order to avoid double counting, we subtract what they have in common, i.e. 
the up-to-C(Q;5) terms in the expansion of Eq.(3.35). Our final result for the 
resummed coefficient function reads: 

V^{mt,mw,Ois{n^),IJi,IJiF) = V%{mumw,Ois{ii^), ii, i^f) 

+ 



as 



as 



(3.36) 



where [f;^]Qs and [f^vjas are respectively the expansion of Eq.(3.35) up to 0{as) 
and the full fixed-order top-decay coefficient function at 0{as), presented in the 
appendix. 

We would like to compare our resummed expression with other similar results 
obtained in heavy quark decay processes [80, 81, 82, 62]. Besides the obvious re- 
placement of a bottom quark with a top in the initial state, our work presents other 
essential differences. We have resummed large coUinear logarithms a5ln(m^/m| 



while Refs. [80, 81, 82, 62] just address the decay of heavy quarks into massless 
quarks. Moreover, this work still differs in a critical issue. Those papers are con- 
cerned with observing the lepton produced by the W decay or the photon in the 
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b Xg'y process, while we wish instead to observe the outgoing b quark. This is 
immediately clear from the choice of the z variable whose 2; — > 1 endpoint leads 
to the Sudakov logarithms. In our case it is the normalized energy fraction of the 
outgoing bottom quark; in [80, 81, 82, 62] it is instead related to the energy of 
either the lepton or the radiated photon. 

The most evident effect of this different perspective is that an additional 
scale, namely the invariant mass of the recoiling hadronic jet, enters the results 
[80, 81, 82, 62], but as we noted above, is absent in our case. An additional 
function (called 'yias) in [80, 81], C{as) in [82], B{as) in [62]) appears in those 
papers. The argument of 0:5 in this function is related to the invariant mass of the 
unobserved final state jet constituted by the outgoing quark and the gluon(s). It 
is worth noting that an identical function, called B [as (Q^(l — z))], also appears 
in the e+e~ [43] and DIS [53] massless coefficient functions, where it is again as- 
sociated with the invariant mass of the unobserved jet. We also observe that to 
order as the coefficient S**-^-* coincides with the corresponding H^^^ of the function 
H [as {ml{l — z)'^)] (see Eq.(2.83)), which resums soft terms in the initial con- 
dition of the perturbative fragmentation function. It will be very interesting to 
compare the functions S{as) and H{as) at higher orders as well. 

One final comment we wish to make is that, as expected, in our final re- 
sult, Eq.(3.22), which accounts for NLL soft resummation in both the coefficient 
function and the initial condition of the perturbative fragmentation function, 

In""^^ N terms do not appear, since they are due to soft and coUinear radiation. 
Both the quarks being heavy, only the former leads to a logarithmic enhance- 
ment. Double logarithms are generated by a mismatch in the lower and upper 
integration limits over the ^[^^(g^)] function in the exponent of the resummation 
expression. In our case both of them have the same functional dependence with 
respect to z, i.e. (1 — zY in Eq.(3.29) and Eq.(2.83) (corresponding to a = 2 in 
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(2.74)). The cancellation of the agin"'*' term can be explicitly seen at order 
as by comparing the large- A?" hmit for the coefficient function, Eq. (3.24), and 
the initial condition (Eq.(A.27)): the In^A^ terms have identical coefficients and 
opposite signs. 

As follows from Eq.(3.22), to obtain the resummed differential distribution in 
X(,-space, one needs to perform the inverse Mellin transformation (2.65). As we 
previously emphasized, we perform that inversion numerically, using the so-called 
Minimal Prescription [63] : the integration contour in the complex A^ plane crosses 
the real line at a value Cmp such that 2 < Cmp < N^. That way all singularities 
of the integrand are to the left of the contour with the exception of the Landau 
pole Nl — exp(l/(2Q!5(/i)6o))- The origin of Nl is in the Landau pole associated 
with the integration over the running strong coupling in the Sudakov exponent 
(3.29), when combined with the relations (2.82), (3.32) and (2.81). 

Next we present the results for the 6-quark energy distribution in top decay 
including the soft-ghion rcsummation in the initial condition of the perturbative 
fragmentation function and the coefficient function. The rest of the settings are 
the same as the ones we discussed in the end of the previous section. 

In Fig. 3-6 we present the xj, distribution according to the approach of per- 
turbative fragmentation, with and without NLL soft-gluon resummation. For the 
scales we have set /x^ = = mt and /xq = A*of = We note that the two 
distributions agree for < 0.8, while for larger Xb values the resummation of large 
terms x;, — > 1 smoothes out the distribution, which exhibits the Sudakov peak. 
Both distributions become negative for — > and x;, — > 1. The negative behav- 
ior at small Xb can be related to the presence of unresummed as In Xb terms in the 
coefficient function. At large Xb, we approach instead the non-perturbative region, 
and resumming leading and next-to-leading logarithms is still not sufficient to cor- 
rectly describe the spectrum for Xb close to 1. In fact, the range of rehabihty of the 
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Fi gure 3-6: 6-quark energy distribution in top decay according to the perturbative fragmenta- 
tion approach, with (sohd hne) and without (dashes) NLL soft-gluon rcsummation. In the inset 
figure, we show the same curves on a logarithmic scale, for xj, > 0.8. We have set fip = 1^ = rrit 
and ijQF = A*o = mi,. 

perturbative calculation has been estimated to be 2:5 < 1 — A/mj ~ 0.95 [43]. It 
is interesting to investigate the dependence of phenomenological distributions on 
the renormaUzation and factorization scales which enter the coefficient function 
{/I and /If) and the initial condition of the perturbative fragmentation function 
{Ho and Hof)- In particular, it is worth comparing the 6-energy spectra with and 
without soft resummation. For the scales fi and fip we consider the values nit/2, 
rrit and 2mt; for fiQp and fiQ the choices are mi,/2, rrib and 2m(,. Figs. 3-7 and 3-8 
show the dependence of the Xb spectrum on the factorization scales /ip and //qf; 
the dependence on the renormalization scales /i and /iq is exhibited in Figs. 3-9 
and 3-10. 

We note that all distributions which include soft-gluon resummation exhibit a 
reduced dependence on the factorization and renormalization scales. 

Fig. 3-7 shows that curves obtained using different values oi iip are almost 



Solid: soft resummed 
— Dashes: non resummed 




CHAPTER 3. B-QUARK FRAGMENTATION IN T-QUARK DECAY 73 




Figure 3-7: 6-quark energy spectrum for different values of the factorization scale fip, with 
(solid) and without (dashes) NLL soft-gluon resummation. The other scales are fixed at fi = mt, 
Mo = MoF = m-b. As in Fig. 1, in the inset figure, we present the same plots for large values of 
Xb, on a logarithmic scale. 




Figure 3-8: As in Fig. 3-7, but for different values of /xqf- The other scales are fixed at 
H = jiF = mt, Ho = mb- 
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Figure 3-9: As in Fig. 3-7, but for different values of the renormalization scale /x. The other 
scales are fixed jxp = mt, = hof = Wfe- 




Figure 3-10: As in Fig. 3-7, but for different values of the renormalization scale /^o- The other 
scales are fixed a.t ^ = fxp = mt, Hof = mb- 
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indistinguishable once soft resummation is included; the unresummed plots exhibit 
a stronger effect of the chosen value ior fip. A similar result also holds for the 
scale /loF'- the dependence of the plots on its actual value for Xf, > 0.8 is small if 
soft logarithms are resummed and quite strong if the prediction is unresummed 
(Fig. 3-8). 

The choice of the value for the renormalization scale n appearing in Eq. (3.24) 
affects only the neighborhood of the Sudakov peak of the resummed predictions, at 
Xb- values very close to one (Fig. 3-9), where, as we have pointed out, our perturba- 
tive approach is anyway unreliable. The effect of the choice of the renormalization 
scale iio on the soft-resummed spectra is slightly larger than the effect of and 
visible at < 1 as well (Fig. 3-10). As for the non-soft-resummed predictions, 
although all dashed curves in Figs. 3-9 and 3-10 seem to converge to the same 
point for x;, — > 1, the overall dependence on /i and /Iq for Xb < 1 is stronger than 
for the resummed predictions. 

As a whole, one can say that the implementation of NLL soft-gluon resumma- 
tion, along with the NLL DGLAP evolution for the perturbative fragmentation 
function, yields a remarkable improvement of our phenomenological results, since 
the reduced dependence on the choice of factorization and renormalization scales 
in the region where the perturbative approach is reliable corresponds to a reduc- 
tion of the theoretical uncertainty. 

3.5 Energy Spectrum of 6-flavored Hadrons in 
Top Decay 

In this section we consider the inclusion of a non-perturbative component in addi- 
tion to the perturbative result, so as to make predictions for observable 6-flavored 
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1 dzdV^ \T^r,nf^B\ 

[xb, mt, mw, mb) = —(2;, rrit, mw, nibjD ^ — , (3.37) 

T JxB z dz \ z J 



Xb = 2~' (3-38) 

1 — w mf 



hadrons (like B mesons) in top decay. At the same time, we also account for the 
inclusion of NLL soft and coUinear resummation. 

We write the normahzed rate for the production of 6-hadrons B as in Eq.(2.64), 
i.e. as a convolution of the rate for the production of b quarks in top decay, given 
by Eq.(3.22) and a non perturbative fragmentation function D'^^[x): 

IdT^ , . 1 dzdV\ 

r dxB i 

where is the B normalized energy fraction: 

1 2pB_Pt 

Pb being the B four- momentum. Since D"'P[x) contains non-perturbative infor- 
mation, it cannot - for the time being - be calculated from first principles in QCD, 
but can only be extracted from data. 

We shall assume a universality property for such a function, i.e. that it does 
not depend on the process but only on the non-perturbative transition it describes, 
and extract it from fits to S-production data collected at LEP in e+e" collisions. 
In particular, we can choose different functional forms for £)"^(x), and tune these 
hadronization models to the available data. We shall consider three models: a 
power law with two tunable parameters: 

0''fe«.ffl = ^(^^^'_^^^ (l-xrx^ (3.39) 

the model of Kartvelishvih et al. [83] : 

L>"P(x; S)^{1 + 5){2 + 5){1- x)x^ , (3.40) 

and the Peterson et al. model [84]: 
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In Eq.(3.39), B{x,y) is the Euler Beta function; in (3.41) A is a normalization 
constant for which an exphcit expression can be found in [45]. 

In order for our extraction procedure to be self-consistent, we employ the 
same underlying perturbative description in both the e'^e" bb process (where 
the non-pcrturbative contribution is fitted) and t — > bW (where it is used). This 
will be ensured by using in both cases NLO, MS coefficient functions, along with 
a fully NLL soft-gluon resummed description, with the large coUinear logarithms 
resummed to NLL accuracy by DGLAP evolution. For the coefficient functions 
in e'^e~ annihilation we shall refer to [85]. 

Fits to data points can be performed either in x^-space, or, as recently advo- 
cated [47], in the conjugated moment space. When fitting in xb space we discard 
data points close to xb — and xb — i and consider ALEPH [86] data in the 
range 0.18 <xb< 0.94. Good-quality data are also available from SLD [87] on 
6-flavored mesons and baryons. Recently, the OPAL [88] and DELPHI [89] Col- 
laborations also published new results, which are fully compatible with the ones 
from ALEPH. The results of our fits are shown in Table 3.1. 

The models in Eqs. (3.39) and (3.40) yield very good fits to the data, while 
the model (3.41) is marginally consistent. We have also fitted the SLD data 
and found qualitatively similar results. However, the values for the parameters 
which best fit the SLD data are different from the ones obtained for ALEPH and 
quoted in Table 3.1. Using the results in Table 3.1 we can give predictions for 
the spectrum of 6-flavored hadrons in top decay. To account for the errors on 
the best-fit parameters, we shall plot bands which correspond to predictions at 
one-standard-deviation confidence level. In Fig. 3-11 we show our predictions for 
the Xb distribution using the three models fitted to the ALEPH data. At one- 
standard-deviation confidence level, the three predictions are different, with the 
Peterson model yielding a distribution which lies quite far from the other two 
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a 


0.51 ±0.15 




13.35 ± 1.46 




2.56/11 


S 


17.76 ±0.62 


x\S)/doi 


10.54/15 


e 


(1.77 ±0.16) X 10-3 


x'(e)/doi 


29.83/15 



Table 3.1: Results of fits to e+e" ^ bb ALEPH data, using matched coefficient function 
and initial condition, with NLL DGLAP evolution and NLL soft-gluon resummation. We set 
A(5) = 200 MeV, hqf = iiq = mb = 5 GeV and /iF=M = \/s = 91.2 GeV. a and /3 are the 
parameters in the power law (3.39), 6 refers to (3.40), e to (3.41). The fits have been performed 
neglecting the correlations between the data points. 

and peaked at larger x^-values. Within two standard deviations, the predictions 
obtained using the models (3.39) and (3.40) are nonetheless in agreement. 

We would like to emphasize that the differences between the various models 
mainly originate from the varying quahty of the fits to e~^e^ data where, as the 

values in Table 3.1 seem to suggest, a given model is sometimes not really able 
to describe the data properly, due to its too restrictive functional form. 

An alternative, and probably better, way of determining and including non- 
perturbative information makes use of moment space perturbative predictions and 
data [47]. The full hadron-level result can be written in A^-space as the product 
of a perturbative and a non-perturbative contribution, = T'j^D^. For each 
value of N one can then extract the corresponding value from e'^e" data, 
with no reference whatsoever to a specific hadronization model, and use it to 
predict the same moment in top decay. The DELPHI Collaboration [89] has 
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Figure 3-11: xb spectrum in top decay, with the hadronization modelled according to a 
power law (solid lines), the Kartvelishvili ct al. (dashes) and the Peterson (dots) model, with 
the relevant parameters fitted to the ALEPH data. The plotted curves arc the edges of bands 
at one-standard-deviation confidence level. NLL soft-gluon resummation is included. We set 
j2p = jjL = rrit and /Uqf = A*o = fnb. 

recently published preliminary results for the moments of 5-meson fragmentation 
in e^e~ collisions up to iV = 5. From these data, and using the moments of the 
e'^e" perturbative contributions [43], one can extract D"^ . The corresponding 
values can then be calculated making use of the results for F^ obtained in [65] 
and given by Eqns.(3.22) and (3.36). CaUing and cr^ the moments for the 
production rate of B mesons (measured) and h quarks (calculated in perturbative 
QCD) in e+e" annihilation, we have cr^ = (t^D^ and hence 

= = ^'n4- ■ (3-42) 

Cat 

Table 3.2 shows a practical implementation of this procedure. Predictions for 
the moments F^ of 5-meson spectra in top decay are given, making use of the 
DELPHI experimental data. Two sets of perturbative results ([A] and [B]) are 
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(x) (x3) (x^) 


e+e data 


0.7153±0.0052 0.5401±0.0064 0.4236zb0.0065 0.3406±0.0064 


e+e- NLL a% [A] 
e+e~ NLL a% [B] 


0.7666 0.6239 0.5246 0.4502 
0.7801 0.6436 0.5479 0.4755 


D7 [A] 
D7 [B] 


0.9331 0.8657 0.8075 0.7566 
0.9169 0.8392 0.7731 0.7163 


t-decay NLL r% [A] 
t-decay NLL L^ [B] 


0.7750 0.6417 0.5498 0.4807 
0.7884 0.6617 0.5737 0.5072 


t-decay [A] 
t-decay [B] 


0.7231 0.5555 0.4440 0.3637 
0.7228 0.5553 0.4435 0.3633 



Table 3.2: Experimental data for the moments cr^ from DELPHI [89], the resummed e+e" 
perturbative calculations for [43], the extracted non-perturbative contribution D^. Using 
the perturbative results F^, a prediction for the physical observable moments is given. Set 
[A]: A(5) = 0.226 GeV and rub = 4.75 GeV, set [B]: A^^) = 0.2 GeV and rub = 5 GeV. The 
experimental error should of course be propagated to the final prediction. 

shown, the first using A^^^ = 0.226 GeV and = 4.75 GeV, the second one using 
our default parameters introduced in Section (3.3). As expected, the perturbative 
calculations and the corresponding non-perturbative components differ, but the 
final predictions for the physical results are to a large extent identical. 



Chapter 4 

Charged— Current Deep Inelastic 
Scattering 

In this Chapter wc present our result for the threshold resummation of the co- 
efficient function for heavy quark production in Charged-Current Deep Inelastic 
Scattering (CC DIS) at low transferred momentum Q (i.e. when Q is of the or- 
der of the mass of the heavy quark) and with next-to-leading logarithmic (NLL) 
accuracy. We then apply our result to the case of charm-quark production in CC 
DIS processes with the masses of the charm quark and the target taken fully into 
account. The presentation in this Chapter is based on the paper [90]. 

4.1 CC DIS: Notation and Overview 

DIS is a process of scattering of a lepton and a hadron (usually nucleon) that 
produces another lepton and final hadronic state X: 

^iniPiJ + N{P) ^ ioutiPi^J + X. (4.1) 

In general, the leptons and lout can be of the same type (neutral-current 
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Figure 4-1: Deep Inelastic Scattering. 

scattering) or be different (charged current scattering). To lowest order in the 
electroweak couphng, the DIS interaction has the structure shown in Fig. (4-1). 
It represents the interaction of a leptonic current j'^ with a hadronic current J'* 
through the exchange of a virtual vector boson. The vector boson is a combination 
of a virtual photon 7 and Z° in the NC case and a virtual charged boson in 
the CC case. 

Let us introduce the following notation: as usual, we will denote by q the 
momentum of the transferred boson (of any type). The momentum is space-like, 
and = —q^ > 0. The momentum of the incoming nucleon is P with P^ — M^, 
M being the mass of the nucleon. We will consider two types of targets: one, 
where is a proton and the other an isoscalar target, where iV is a "half-sum" of 
a proton and a neutron. In the latter case, the partonic content of is determined 
by assuming isospin invariance (2.3). We will not distinguish between the masses 
of the two target types. Numerically, in our applications we set M = 1 GeV. A 
very important parameter is the Bjorken variable: 

which is directly measured in the experiment. That variable has the nice property 
that in the parton model it coincides with the momentum fraction carried by the 
initial state parton, as was discussed in Section (2.1.2). Another common variable 
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is: 

- = (4.3) 

as well as the inelasticity y: 

Prom the discussion following Eq.(4.1) it becomes clear that the cross section 
for DIS process factorizes into leptonic and hadronic parts [91]: 

da ~ L^'W^,, (4.5) 

where: 

V = ^'Zi^in\jl\iout){iout\mn), (4-6) 

is a tensor containing the leptonic current and 

W^.^j^Y.i'^7r)'^\P + Q-Px){P\MPx){Px\Jl\P), (4.7) 

is the tensor containing all the hadronic contributions. 

It is obvious that unlike L^,^, its hadronic counterpart W^i, cannot be simply 

calculated, since the latter contains non-perturbativc information incorporated 
in the hadronic states present in Eq.(4.7). It is customary to parameterize the 
hadronic tensor in the following way: 

W,. = -g,.W, + - ^'-^^^^f^W, + ..., (4.8) 

where the functions 14^1,2,3 are scalar form-factors, and ". . ." stands for other 
contributions that do not contribute in the limit when the lepton masses ^ are 
neglected. It is however more common to work in terms of the re-scaled form- 
factors: 



^For more details on the treatment of lepton masses see [91] and especially [92]. 
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F3 = ^W,. (4.9) 

One can then write the cross-section in Eq.(4.5) for the case of neutrino-nucleon 
scattering (which is what we need for our apphcations) in the following form: 



dxdy 7r(l + gVM^)2 



y'xF, + [l-[l + -^JyjF2±y[l-^]xF,, 



Mx\ \ „ . / y 



(4.10) 

where E^, is the neutrino energy in the nucleon rest frame and y is given by (4.4). 
Similar results exist for the NC case. The functions ^1,2,3 are known as structure 
functions. The knowledge of these is equivalent to the knowledge of the DIS 
cross-section. Since the structure functions are observables they are free from 
ambiguity. In the following we turn our attention to their evaluation. 

Following [91], and in accordance with our discussion in Section (2.2) of the 
factorization theorem, one can write the DIS cross-section as a convolution of a 
parton distribution and a partonic cross-section: 

We prefer to work in terms of the tensor W^^i, rather than da since the former 
already contains all the relevant information for the hadronic part of the reaction. 

The partonic tensor Wni, is defined in complete analogy with W^,^ (cf. Eq.(4.7)), 
but with the hadron target replaced by a partonic one. A parametrization similar 
to Eq.(4.8) holds for w^jy as well. One just needs to replace in (4.8) the hadron 
momentum P with the momentum p of the incoming parton a, and the factors 
containing H^2,3 are normalized with respect to Q'^ instead of M^. The dependence 
of the partonic tensor W/^i, on the convolution parameter ^ (which is the energy 
fraction carried by the parton a) is through the relation p+ = between the 



W^.{q,P) jfa/N{i,lAK^{q,P,ocs{lA)- (4.11) 
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plus light-cone components of the partonic and hadron momenta (see also the 
discussion following Eq.(2.22)). We will discuss that point in more detail below. 

As follows from our general discussion in Section (2.2) and Section (2.4.1), the 
partonic tensor w^i, can be calculated in perturbation theory to any order in the 
strong coupling. As it was explained there, in the case of massless initial-state 
partons ^ the tensor w^i, contains a coUinear divergence. That divergence must 
be subtracted (i.e. effectively absorbed into the parton distribution functions), 
and one usually does that in the MS scheme. This is also the subtraction scheme 
that we use in our applications. The partonic tensor w^i/ depends on the parton 
momentum p through the partonic analogue x of the Bjorken variable, x is defined 
by Eq.(4.2) but with the momentum P replaced by p. 

In fact it is easy to see from the kinematics of the partonic scattering process 
for production at NLO of quark with mass m 

qi{p)W*{q) ^ q2{p2){+X), (4.12) 

that the partonic cross-section can be evaluated in terms of the following variable: 

X 2p-q ' - - ^ ^ 

In Eq.(4.13) we have also introduced the following parameter: 

Although the tensors W^i, and w^i, have the same tensor structure, they are 
not simply proportional to each other because of their dependence on different 
momenta (P and p respectively). Therefore the relation between the invariant 
hadronic and partonic form-factors Wi and Wi depends on the relation between 



^In our discussions we will only consider the case of massless incoming partons. Results for 
the case > can be found in [91]; see also [41]. 
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the momenta P and p. One can parameterize the momentum p of the parton a 
in the following way [91]: 

= cpP'' + Cgq^", (4.15) 
where the coefficients cp and Cq read: 

In Eq.(4.16) ^ is the momentum fraction carried by the incoming parton and rj 
is the Nachtmann variable [93], that corrects the Bjorken variable x for non-zero 
target mass M: 



1 1 / 1 M2 1 1 M2 ^ ^ 



It is now obvious that the whole dependence of the partonic cross-section Wfj_,^{z) on 
the momentum fraction ^ (that is also the convolution variable) enters implicitly 
through the relation (4.15). One can make this dependence explicit by combining 
(4.15) with the other relations above. After some algebra one can show that: 

^ = (4.18) 

where, as is common in the DIS literature, we have introduced the following 

parameter: 

X-j- (4.19) 

The imphcation of the relation (4.18) is clear: one can calculate the partonic tensor 
in perturbation theory as a function of the variable z (as defined in Eq.(4.13)), 
and then use the relation (4.18) to show that indeed the product of the two 
terms in the integrand of Eq.(4.11) is of the general form (2.20). From simple 
kinematical considerations (i.e. using that Spartonic > "^2) one can also show that 
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X < ^ < 1. Clearly, that way the hadronic tensor PV^j, in essence becomes a 
function of the variable x- In the limit of vanishing quark and target masses, that 
variable coincides with the usual Bjorken variable x: x — > x. 

It will be, however, more convenient for us to work in terms of the structure 
functions because the w^,^ depends on ^ not only through the arguments of the 
scalar form-factors Wi_2,3, but also through the partonic momenta defining the 
tensor structure of w^i, (cf. Eq.(4.8)). In the following, in order to relate the 
hadron level structure functions with their parton level counterparts, we follow 
the notation in [94] and [92] . First one defines the "theoretical" structure functions 
J-'i which are related to the Fj's of Eq.(4.9) and (4.10) via the following relations: 

F^ = (4.20) 

2x 

F, = —T, (4.21) 

= -^3, (4.22) 
P 

with: 



The introduction of the structure functions jFj is convenient since they can be 
straightforwardly expressed as a convolution of parton distribution functions and 
MS coefficient functions: 

+ Cf ^^,ij\^il,X^ g{^,^il) , ^ = 1,2,3. (4.24) 

In Eq.(4.24), we have exphcitly shown the relevant arguments of the structure 
functions, although it is clear that the dependence on x comes through the variable 
X in the following way: 
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That relation follows from Eqns.(4.17) and (4.19). Bjorken x is constrained to be 
in the range: 

As usual, in Eq.(4.24) /j, and hf are the renormalization and factorization scales; 

and are correspondingly the coefficient functions for the quark-scattering 
subprocess 

qi{p)W*{q)^q2{p2){9{Pg)), (4.27) 
and for the gluon-fusion sub-process 

9{Pg)W*{q)^UpMp2). (4.28) 

In Eq.(4.24), by qi{x^^^p) and g{x,^\) we have denoted the parton distribution 
functions of the initial-state light quark and the gluon. For consistency, one has 
to use MS parton distribution functions. 

4.2 Behavior of the Coefficient Function in the 
Soft Limit 

In our applications, we would like to consider 0{as) corrections to CC DIS. The 
coefficient functions can be evaluated perturbatively and to 0{as) they read: 

C!{z,i,\\) = Sil-z) + ^4^Hf{z,i,l,X) (4.29) 



27r 

The exphcit expressions for the functions i^i'f 3; with the coUinear divergences 



Cfiz,^,',^^l,X) = ^i^Hfiz,f^%,X). (4.30) 



subtracted in the MS factorization scheme, have been obtained in [95] and [94]. 
For our applications, we use the results for iff f 3 reported in 
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As we are inclusive with respect to the final-state heavy quark q2 the quark- 
scattering coefficient functions are free from logarithms ln(m^/(5^). However, the 
functions Hf contain terms which behave like ~ 1/(1 — 2;)+ or [ln(l — z) / {1 — z)]+ 
[94] and are therefore enhanced once the quark energy fraction approaches one. 
The z ^ 1 limit corresponds to soft-gluon emission. In Mellin moment space, such 
a behavior corresponds to contributions ~ InA'^ or ~ InA'^^ in the limit N ^ oo; 
see the appendix for more details. 

The gluon-initiated coefficient functions Hf are not enhanced in the limit z — > 
1, since the splitting g ^ qq is not soft divergent. Instead, they contain terms 
~ ln'^(l — z) of collinear origin but such terms are suppressed in Mellin space by 
an inverse power of N and are therefore neglected at our NLL level of accuracy. 
The gluon- fusion coefficients Hf also contain coUinear logs ln(m^/(5^) originating 
in collinear gluon splitting. The treatment of those terms is related to the choice 
of the parton distribution. We will return to that point later in this Chapter. To 
summarize, in the rest of our discussions we will not be concerned with the soft 
limit of the gluon-initiated sub-process. 

We would like to make further comments on the light quark initiated coefficient 
functions at large z. Omitting terms that are not enhanced in the soft limit for 
any value of the mass ratio m'^/Q'^, one has: 



Hf{z ^ 1, f^l, A) = H^"f\z, f^l, A) (4.31) 



where: 



1 1-z 

+ -7 TZ ^1 + 



In K 



. (4.32) 



4 V(1-A2)V+ (1-^) 
One immediately sees that the behavior of the coefficient function (4.32) at 
large z strongly depends on the value of the ratio m/Q. To illustrate that point. 
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we consider the two extreme regimes: the small-m/Q hmit, i.e. A ~ 1, and the 
large-m/(5 hmit, i.e. A <^ 1. 

In the first case, setting A = 1 in (4.32) one recovers the large- A?" hmit of the 
MS massless coefficient function reported in [58]: 

as{n'^)CF [1,2 



r<soft I 

I 



=1 = 1 + 



TT 



-WN^ 



4 



\nN\ 



(4.33) 



where 7^ = 0.577 ... is the Euler constant. 

In the second case, for A — > and 2; — > 1, one has Az = A in (4.32). As a 

result, the term ~ [(1 — z)/{\ — \zY]j^ is regular in the soft limit. The following 
expression holds: 



(jsoft 



N |A<1 



a5(A*^)CF 



TT 



27E + 1 - In 



In AT 



(4.34) 



For later convenience we have introduced the following scale: 



= f 1 + 



(4.35) 



In evaluating the Mellin transforms, we have made use of the results presented in 
the appendix. 

It is obvious that the coefficient functions are enhanced in the limit z ^ 1 
and therefore must be resummed. Let us note, however, the very different behavior 
of the NLO result in the soft hmit for small and large values of the ratio m/Q (see 
Eqns.(4.33) and (4.34)). That suggests that the resummed coefficient function 
may exhibit strong dependence on the value of m/Q. In the next section we 
will present our results for the resummation of the coefficient function with NLL 
accuracy for low transferred momentum Q"^ ~ m^, and will discuss its relation to 
the well known massless limit vr? jQ^ — > 0. 
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4.3 Soft-gluon Resummation for the Quark-initiated 
Coefficient Function 

One can perform the soft-gluon (threshold) resummation with NLL accuracy for 
the coefScient function for the process (4.27) by applying the general methods 
developed in [53] that we presented in detail in Section (2.4.3). To that end, one 
first needs to collect the contributions from soft-gluon emission. That can be done 
using the eikonal techniques discussed also in Section (2.4.3). Second, one needs 
to identify and - if it is present - to collect the contribution from pure coUinear 
radiation. 

Let us start with the calculation of the soft contributions in the eikonal ap- 
proximation. We consider the inclusive hard scattering process 

qi{p)W*{q) 52(^2) + soft gluons, (4.36) 

where the momenta p, q and p2 are hard and the momenta of the additional gluons 
are all soft. In cm. co-ordinates we have: 

2p0 ^ ^ . (4.37) 

^(l-z)Q2 + ^2 

In the eikonal approximation the momentum of the outgoing hard quark is taken 
[53] to he p = p + q, with p^ given by: 

f c:^{l-z){Q^ +m^)+m^. (4.38) 

Next we calculate the contribution from real gluon emission to the eikonal cross- 
section at order 0{as) and add the virtual contribution using unitarity, as in 
(2.73). The eikonal current (2.72) is constructed out of the hard momenta p and 
p, and its square is given in [53] . We then express the result as an integral over the 
variables z and = {p + Pgf{l - z) with < A;^ < {2p°y{l - zf and < z < 1. 
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Finally, the exponentiated result for the coefficient function for the process (4.36) 
in the soft limit and in the MS-subtraction scheme reads: 

In Ar, = dz '\\~^'^ | fj^"" ^A [asi^)] + S [as {kijj] | . (4.39) 

The function A is defined in (2.77) and the function 5" in (3.30). To obtain the 
latter we used the relation ^ (2.75). The derivation of (4.39) proceeds in complete 
analogy with our discussions in Section (2.4.3) and Section (3.4), so we do not 
present each step in full detail here. In the present case, however, where the 
kinematics of the process we consider allows the ratio of the two mass parameters 
m? /Q"^ to vary widely, wc encounter another complication that we discuss next. 

Ideally, one would like to be able to describe the soft limit z — > 1 of the coef- 
ficient function (4.29) with a single expression that is valid for any kinematically 
allowed value of the ratio m^/Q^. It principle that seems possible because the 
zero-mass limit — > of the coefficient function is regular for any value of 
the kinematical variable z. This follows from the fact that the inclusive quark- 
initiated coefficient function is free from mass singularities from the final state. 
However, as can be explicitly seen from the fixed order result we discussed in the 
previous section, the limit 2; — > 1 of the coefficient function is strongly affected by 
the value of rr? jQ"^. The reason is obvious: if one wants to be able to cover all 
values of ni^ /Q"^, effectively one has to deal with a two variable problem. As such 
one can consider z and \z or, in the conjugated Mellin space, N and (1 — A)A^ 
(recall the definition of A in Eq.(4.14)). 

Let us demonstrate our point in a simple example. Consider a term of the form 
ln[(l — A)iV]. We are interested in its behavior at large for any A < 1. Although 
such a term is formally divergent when N ^ 00 for any A < 1, it is clear that 
for values of A sufficiently close to one, such a term cannot be considered large; 
^ After Eq.(4.40) has been used as well; see below. 
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in particular it cannot be considered on the same footing with pure In(A^) terms. 
Another example is the scale /c^ax- behavior strongly depends on the relative 
magnitude of with respect to {1 — z)Q'^. For small enough mass, that scale is 
I — z) while for sufficiently large mass it is ~ (1 — z)^: 



g2(i - z) when = 0, 

^ ^ (4.40) 

M^il-zf whenm^^g^ 



2 

max 



with Ai defined in Eq.(4.35). That difference in the behavior of the scale k 
influences significantly the soft limit of the coefficient function; we will see that 
the coefficient of the leading logarithmic term in the perturbative result (4.33) 
and (4.34) is in fact proportional to the value of the power of (1 — 2;) in /Cj^ax- Yet 
another illustration of the interplay between the large N limit and the value of the 
mass m (relative to Q) is the result for the factor K (4.47), (4.58); the derivation 
of this factor is presented in some detail in the appendix. 

In order to avoid the above mentioned subtleties and to be able to extract 
in a straightforward way the singular z dependence of the resummed coefficient 
function (4.39), in what follows we will consider separately the following two cases 
(each case effectively involves a single scale): the case of sufficiently large value 
of the ratio m^/Q^ that we will refer to as "massive", and the "massless" case 
VP? IQ^ — 0. We will present separately results for both cases. We will treat 
these two results as complementing each other; we will apply the "massive" result 
for cases where w? /Q"^ ^1/5 while we can certainly apply the massless result for 
nn?' IQ^ <1/1Q. Supposedly in the intermediate region neither of the results is 
correct although the difference between the two results in that range is not large, 
i.e. the correct results also should not differ much. 

Finally, in order to collect all contributions one must also account for pos- 
sible coUinear radiation from the final state. As we will see in a moment, that 
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contribution is also affected by tlie value of /Q"^. 



4.3.1 The case w? ~ 

If the final-state quark is massive, the virtuality of the final state in the large-^ 
limit is given in Eq.(4.38) and is never small provided that m/Q is sufficiently 
large. As a result, there is no additional coUinear enhancement for heavy quark 
production. Therefore in the case w?/Q'^ ~ 1, by combining Eqns.(4.39) and 
(4.40), we get the following result for the quark-initiated coefficient function in 
the limit z ^ 1: 

InA^ = / dz— \ l —Aas{k^) 

Jo 1 — z [J^ij, L -I 

+ s[as(M\l-zf)]^. (4.41) 

One can evaluate Eq.(4.41) with NLL accuracy with the help of Eq.(2.82). The 
result has the general form (2.80): 

A^(a5(/^'),/^',/^'F,^',Q') =exp[ln7V^W(L) + ^(2)(^^^2^^2^)j , (4-42) 

where similarly to Eq.(2.81) we have introduced: 

L = hQas{ii^)\nN . 



The functions g^^^ and g^"^^ are given by 

^(1) 



+ 



27rboL 
27r6o 



[2L + (1 - 2L) ln(l - 2L)] 



(4.43) 



(4.44) 



In 



27i 



ln(l - 2L) 



+ 



4L + 2 ln(l - 2L) + \n^{l - 2L) 
-— [2L + ln(l - 2L)] — - + In ^ 
ln(l-2L). 



27r6i 

^(1) 



2nbo 



(4.45) 
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In Eq.(4.42) the term In Ng^^^L) accounts for the resummation of leading loga- 
rithms (LL) In'*"'"^ N in the Sudakov exponent, while the function gr(^) (L, /i^, fip) 
resums NLL terms ccg In" N. 

Next, we present our final Sudakov- resummed coefficient function with matched 
NLO and NLL contributions: 



1 + 



X exp 



— Ki{iJ,p,m ,Q ) 

ln7V^W(A) + ^(^)(A,A^,/Xf) 



(4.46) 



Furthermore, in analogy with the top decay case we considered in the previous 
Chapter and following [43], we include in our final Sudakov- resummed coefficient 
function also the terms i^(/Xp, m^, Q^) of the NLO coefficient function (4.29) which 
are constant in the limit N ^ oo: 



+ A- 2U2 - 



+ 2(7^-l)(7s + 2). 



2A 



(4.47) 



The evaluation of the factor (4.47) is detailed in the appendix. Following 
have defined the quantities Ai, which satisfy the following relations: 



= A3 = 0, A2 = -(1 - A) ln(l - A). 
A 



we 



(4.48) 



We now match the resummed coefficient function to the exact first-order result: 
we add the resummed result to the exact coefficient function and, in order to avoid 
double counting, we subtract what they have in common. The final expression for 
the resummed coefficient function reads: 

^-^(/.^/.^,m^g^) = i/^,(/x^/x^,m^Q^) 

- HN,i{lJ>^,lJ>F,m^,Q^) 

I J as 

+ C%Aii\iil,m\Q^) (4.49) 
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where ^ is the MeUin transform of the perturbative result for the functions Cf 
introduced in (4.29). 

Before closing the discussion of the threshold resummation in the massive 
case we would hke to discuss a relation between the result Eq.(4.41) and the 
corresponding results for other processes. 

One can relate our result Eq.(4.41) to the one for heavy quark decay. An exam- 
ple is the top quark decay Eq.(3.3) we discussed in the previous Chapter. One can 
easily see that the expression Eq.(4.41) can be obtained from the corresponding 
result in top decay - Eq.(3.29) - after the identification: 

mtop rn ; -Q^ . (4.50) 

The relation (4.50) is not accidental. It can be understood in the following 
way: consider the coefficient functions for the following two processes: 

Q^W^q, (4.51) 
g + (4.52) 

By Q and q we denote heavy and light quark, respectively, and the states in 
Eqns.(4.51), (4.52) are the initial/final states in those processes ^. Ignoring details 
like spins and charges, the processes (4.51) and (4.52) are related ^ by time- 
inversion which, in turn, implies (4.50). Let us now generalize those processes 
and consider more particles (denoted generally by X) in the final states: 

Q ^ + g + X, (4.53) 

q + W ^ Q + X. (4.54) 

^Clearly, that is unphysical in view of the Kinoshita-Lee-Nauenberg Theorem discussed in 

Section (2.4). However, one can think of Eqns.(4.51) and (4.52) as the Born approximation to 

the corresponding physical processes. 

^Plus inversion of the sign of the mass of the Vl^-boson since the latter is space-like in (4.52) 

and time-like (and on-shell) in (4.51). 
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Of interest to us are processes where X stands for a collection of soft gluons. In 
that case, up to contributions that are not singular in the limit — > oo, the 
processes (4.53) and (4.54) can also be related by Eq.(4.50). To see this, one 
needs to recall that the processes (4.53) and (4.54) can be described with the help 
of the eikonal cross-section discussed in Section (2.4.3). In this approximation 
the "back reaction" from soft gluon radiation on the hard radiating particles is 
neglected. To derive the soft coefficient function in the eikonal approximation, 
we integrate Eq.(2.68) over all of the radiated gluons. The resulting expression 
therefore depends only on the configuration of the hard scattering particles that 
in turn is encoded in Eqns.(4.51) and (4.52). 

An explicit example is the one-loop evaluation of the eikonal amplitudes in our 
processes of interest: the results depend only on the square of the eikonal current 
for which (4.50) is manifest. 

Of course the above considerations are complete only because there is no ad- 
ditional contribution due to collinear radiation in any of the processes (4.53) and 
(4.54). We would also like to note that the relation (4.50) is not restricted only to 
the NLL approximation. For that reason we expect that the function S appearing 
in both Eq.(4.41) and Eq.(3.29) is the same even beyond the NLL level. 

Let us point out that soft resummation for the coefficient function in DIS has 
been first considered in [53, 58] and lately implemented in [61], but in the approx- 
imation in which all participating quarks are massless. Soft-gluon resummation 
for heavy quark production in DIS processes has been investigated in [96] and 
[97] , where the authors have considered neutral current interactions. 
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4.3.2 The case < 

In the limit m^/Q'^ ^ we recover the well known soft limit of the massless 
coefficient function (4.29) derived in [53]. As follows from Eq.(4.38), in this case 
the virtuahty of the final-state jet vanishes hke (1 — z)Q'^. Therefore the NLL 
coefficient function receives an additional contribution due to purely coUinear 
radiation. As shown in [53], one can account for such a contribution by solving a 
modified evolution equation for the jet distribution associated with the final state. 
As a result, an additional term will contribute to the function S{as) that enters 
the Sudakov exponent in Eq.(4.39). Taking also into account the zero-mass limit 
in Eq.(4.40), we get the well known massless result for the quark initiated MS 
coefficient function [58]: 

fi f rQH^-^) dp r i 1 r / \il 

In AnU/q^o = dz—^—— I j ^ —A [as{k^)\ + -B [ag - z))\ \ . 

(4.55) 

The function B{as) generalizes the function S from Eq.(4.39) and also includes the 
contribution from coUinear radiation from the final state. It also can be expanded 
in powers oi ag'- 

5(«5) = f:(— (4.56) 

n=l ^ ^ ^ 

and, to NLL level, one needs only the first term of the expansion: 

= J-Cp. (4.57) 

The explicit evaluation of Eq.(4.55) to NLL level can be found in [61], where 
functions analogous to our g'^^'^ and g'^'^^ are reported. Similarly to Eq.(4.49) in the 
massive case, we also match the resummed result (4.55) to the exact coefficient 
function (4.29), with functions now evaluated in the limit m/Q — > 0, i.e. 
A — > 1. The constant factors Ki, analogous to the one in Eq.(4.47), read in the 
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massless approximation (see the appendix for details): 

K^ii^l, m^ Q') U/Q^o ={l- 2^e) In + ^| + - !^ - | (4.58) 

for i=l, 2, 3. 

Finally, one can check that the 0{as) expansion of the resummed results 
correctly reproduces the large- limits of the coefficient function in both the 
massive and massless cases. 

Before we move on to the phenomenological applications, let us compare our 
results on soft-gluon resummation for heavy quark decay Eq.(3.29) and for heavy- 
quark production in CC DIS Eq.(4.41) to the general results for soft-gluon resum- 
mation discussed in [21]. Both our results have the general form: 



lnA5^^ = lnA^ + lnJ^, (4.59) 



where the factor: 



is due to soft and coUinear radiation from the final quark for the case of top-decay 
(similarly to the e'^e" quark production in [43]) and from the initial quark in CC 
DIS. The scale appearing there is the corresponding characteristic hard scale. 
It equals mj{l — rn^/mlY' for the case of top-decay and m^(l + Q"^ /in^Y fo'^ the 
case of massive quark production in CC DIS. The second term in Eq.(4.59) 

InJ^ = j\z^-Y—^S[as{Q\l-zf)\, (4.61) 

describes soft-gluon radiation that is not coUinearly enhanced, and the scale is 
the same as the one in (4.60). Our result (4.61) however differs in an obvious way 
from the similar factors appearing in [21, 43]: in both of our cases (top decay and 
CC DIS) we have an observed massless quark and an unobserved massive quark. 
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For that reason in Eq.(4.61) we have the function S instead of the function 1/2B, 
where the latter embodies coUinear radiation from a massless quark. Also, the 
integral over the (eikonal) anomalous dimension that appears in those papers is 
also absent in our cases, since for us the two limits of that integral coincide, i.e. 
in the massive case we have: 



where a — (1,2) and J — {1/2B, S) in the massless and massive cases respectively. 
Again, Q stands for the relevant hard scale for the process under consideration. 

4.4 Phenomenological Results for Charm Quark 
Production 

In this section we apply our results for soft-gluon (threshold) resummation to the 
study of charm quark production in charged current (CC) events. Measurements 
of the charm quark structure functions in the CC regime are in fact important in 
order to probe the density of strange quarks and gluons in the proton. 

Charged-current Deep Inelastic Scattering processes U/j^N nX were recently 
studied at the NuTeV experiment at Fermilab, where neutrino and antineutrino 
beams coUide with an iron fixed target (see, e.g., Ref. [98] for the updated re- 
sults). In particular, production of two oppositely-charged muons at NuTeV [99] 
is mainly associated with a CC event with a charm quark in the final state. The 
first generation of the HI and ZEUS experiments on electron-proton DIS at the 
HERA coUider at DESY (HERA I) did detect charged-current events ep 
[100, 101], but it did not have sufficiently-high statistics to reconstruct heavy 



InJ^ 




(4.62) 
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quarks in the CC regime. However, such measurements are foreseen in the cur- 
rent improved-luminosity run (HERA H) and at the possible future generation 
of HERA HI experiments [102], which should be able to investigate the proton 
structure functions also at low values. 

The structure functions are given as a convolution of MS coefficient functions 
and parton distribution functions. Since the resummed coefficient function is 
given in A'"-space, in principle one would like to use parton distribution functions 
in Mellin space as well, in order to get the resummed structure function in moment 
space and finally invert it numerically from N- to x-space. However, all modern 
sets of parton distribution functions [29, 103, 26] are given numerically in the form 
of a grid in the {x, Q^) space. 

In order to overcome this problem, we follow the method proposed in Ref. [104] 
in the context of joint resummation, which allows one to use x-space parton dis- 
tributions, even when performing resummed calculations in Mellin space. 

This method consists of rewriting the integral of the inverse Mellin transform 
of the resummed structure functions as follows: 



where Qn^IJ^f) parton distribution of the initial-state quark gi, as if it were 

available in space, and F^v is the integration contour in the complex plane, 
chosen according to the Minimal Prescription [63]. 

One can integrate out the term with the parton distribution function observing 
that the following relation holds: 





(4.63) 



2m 



1 
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Due to Eq.(4.64), one no longer needs the iV-space parton distributions, but one 
should just differentiate the x-space ones, which is numerically doable. 

As one has the analytical expression for if^^^ in A^-space (4.49), the following 
inverse Mellin transform is straightforward: 

= 2^? L diVr" (4.66) 

which, thanks to the suppressing factor 1/[N — lY in the integrand function, turns 
out to be smooth for ^ — > 1. 

The resummed structure function will be finally expressed as the following 
convolution 

After having clarified how we are dealing with the parton distribution func- 
tions, we are able to present results for soft-resummed structure functions. We 
shall consider charm quark production, i.e. ^2 = c in Eq.(4.27), since processes 
with charm quarks in the final state play a role for structure function measure- 
ments at NuTeV and in HERA experiments. 

As for the choice of the parton distribution set, in principle, once data on 
heavy quark production CC DIS were to be available, one should use the CC NLL 
coefficient function when performing the parton distribution global fits and get 
NLL parton densities as well. For the time being, we can just use one of the most- 
updated NLO sets and convolute it with the fixed-order or resummed coefficient 
function. 

We shall present results based on the new generation of CTEQ NLO MS parton 
distribution functions [29], the so-called CTEQ6M set, but similar results can be 
obtained using, e.g., the MRST [103] or the GRV [26] sets. 

The elementary scattering processes which yield the production of charm 
quarks in CC DIS are dW* c and sW* c. For e+(e~)p scattering at HERA, 
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our parton distribution function qi{i,Q^) in Eq.(4.24) will hence be: 

gi(e,Q')|HERA = \VcdM^,Q') + |v;,ps(c,g2), (4.67) 

where Vcd and Vcs are the relevant Cabibbo-Kobayashi-Maskawa matrix elements. 
For neutrino scattering on an isoscalar target, that will have to be modified in 
order to account for the possibility of an interaction with a neutron as well: 

where d{^, Q^), u{^, Q'^) and s{^, Q'^) are still the proton parton distribution func- 
tions and we have applied the isospin symmetry (2.3). 

In order to be consistent with the use of the CTEQ parton distribution func- 
tions, we shall use for the QCD scale in the MS renormahzation scheme the values 
A4 = 326 MeV and A5 = 226 MeV, for four and five active fiavors respectively. 
This corresponds to as{mz) = 0.118. The charm and bottom quark masses have 
been set to nic = 1.3 GeV and rrih = 4.5 GeV, as done in [29]. 

Since mass effects are important for large values of the ratio rric/Q, we would 
like to have Q as small as possible to be able to apply our massive resummation. 
The NuTeV experiment is indeed able to measure structure functions at small Q"^ 
values [98, 99], as they are able to detect the final-state muon produced in the 
p^N fiX processes. In our phenomenological analysis we shall consider charm 
production at = 2 GeV^ and Q'^ — 5 GeV^, which are values reachable at 
NuTeV and such that ml/Q'^ is relatively large. 

The detection of CC events at HERA is more problematic, due to backgrounds 
and to the presence of a neutrino, instead of a charged lepton, in the final state. 
The current HERA 11 experiments make it possible to detect heavy quarks in CC 
events, but the values for such events are still supposed to be much larger 
than the charm quark squared mass, typically > 100 GeV^. We shall present 
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Figure 4-2: Results for F2{x, Q^) for charm quark production in neutrino scattering at = 
2 GeV^ with (sohd) and without (dashed) soft resummation in the coefficient function. We have 
set = 1^ = Q- 

nonetheless results for charm quark production at HERA, but in this case we 
shall have to use the massless result for the resummed MS coefficient function, i.e. 
Eq.(4.55) and the formulas reported in [61]. We shall consider the typical HERA 
values = 300 GeV^ and 1000 GeV^. 

We shall present results for the structure function ^2^(0;, Q^), but we can al- 
ready anticipate that the effect of the resummation is approximately the same on 
all three structure functions , F| and Fg and on the single-differential cross 
section da/dx, obtained from Eq.(4.10), after integrating over y. 

For the sake of comparison with the experiments, we shall plot the structure 
function F| in terms of Bjorken x, which is the quantity which is measured. For 
most of our plots, we shall set the renormalization and factorization scales equal 
to Q, i.e. /If — fj, — Q. Afterwards, we shall also investigate the dependence of 
our results on the choice of these scales. 




CHAPTER 4. CHARGED-CURRENT DEEP INELASTIC SCATTERING 105 



u 

CM 

10-4 
10-^ 

0.0 0.2 0.4 0.6 0.8 

X 

Figure 4-3: As in Fig. 4-2, but for = 5 GeV^ 

In Figs. (4-2) and (4-3) we show F|(a;, Q^) in neutrino scattering at — 
2 GeV^ and 5 GeV^ and compare the results obtained using fixed-order and soft- 
resummed MS coefficient functions. We have set the target mass to M = 1 GeV, 
which is a characteristic nucleon mass. Wc observe a relevant effect of the imple- 
mentation of soft resummation: the two predictions agree up to a; ~ 0.2 — 0.3; 
afterwards one can see an enhancement of the structure function due to the resum- 
mation. At very large x resummed effects are indeed remarkable: for Q'^ — 2 GeV^ 
one has an enhancement of a factor of 2 at x = 0.5 and a factor of 5 at x = 0.6. 
For = 5 GeV^ one gets a factor of 2 at a; = 0.7 and 5 at a; = 0.8. 

In Figs. (4-4) and (4-5) we present results for F|(a;, Q'^), but for charm produc- 
tion at HERA, in particular for positron-proton scattering at = 300 and 1000 
GeV^. In this case, since rric/Q <^ 1, we use the massless result (4.55) for the 
resummed coefficient function. We observe that the impact of the resummation 
is smaller than in the case of low values. This is a reasonable result: in fact, 
leading and next-to-leading logarithms in the Sudakov exponent are weighted by 
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0.0 0.2 0.4 0.6 0.8 1.0 



Figure 4-4: Results for (^i Q^) ^'^^ charm quark production in positron-proton scattering at 
HERA, for = 2 GcV^ We have set /j-f = IJ^ = Q ■ In the inset figure, we show the same plots 
at large x and on a linear scale. 




Figure 4-5: As in Fig. 4-4, but for = iqoo GeV^. 
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Figure 4-6: Dependence of F2 on the factorization scale for neutrino scattering at NuTeV. 
Solid lines include soft resummation in the coefficient function, dashed lines are fixed-order 
predictions 

powers of Q;s(/.i^). As, e.g., 0:5(2 GeV^) ~ 3 as(300 GeV^), resummed effects are 
clearly more important when is small. Moreover, the larger is, the larger 
the values of x are at which one is sensitive to Sudakov effects. 

We note in Figs. (4-4) and (4-5) that, at about x > 0.6, the fixed-order and 
resummed predictions start to be distinguishable. We estimate the overall impact 
of large- ,T resummation on F^{x, Q^) at = 300 GeV^ and = 1000 GeV^ to 
be between 10% and 20%. 

In Figs. (4-6) and (4-7) we investigate the dependence on the factorization 
and renormalization scales. We still keep /i — /ip, which we allow to assume 
the values Q^, 2Q'^ and 4Q^. We consider values of 2 and 300 GeV^ for 
the experimental environments of NuTeV and HERA respectively. We see that 
the curves which implement soft resummation in the coefficient function show 
a weaker dependence on the chosen value for the factorization/renormalization 
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Figure 4-7: As in Fig. 4-6, but at HERA for = 300 GeV^. In the inset figure, we plot the 
same curves at large x and on a linear scale. 

scales. In Fig. (4-6) one can see that one still has a visible effect of the value of 
such scales, but the overall dependence on np and fi of the resummed prediction 
is smaller than for the fixed-order. The plots at large exhibit in general a weak 
dependence on /ip and /i even at NLO, as shown in Fig. (4-7). However, while 
the NLO structure function still presents a residual dependence on the scales, 
the resummed predictions obtained with three different values oi fip are basically 
indistinguishable. A smaller dependence on such scales implies a reduction of the 
theoretical uncertainty of the prediction and is therefore a remarkable effect of 
the implementation of soft gluon resummation. 

We finally would like to compare the impact that mass effects and soft-gluon 
resummation have on the charm structure functions. To achieve this goal, we 
plot in Fig. (4-8) the theoretical structure function ^2 (X' defined in Eq.(4.24) 
as a function of the variable x (see Eq.(4.19)), for neutrino scattering at — 
2 GeV^. We compare fixed-order massive (dashed line), fixed-order massless (dots) 
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Figure 4-8: Comparison of massive (dashes) and massless (dots) fixed-order ealculations with 
the resummed massive result (solid) for = 2 GeV^ and fi = jip = Q- Plotted is the theoretical 
structure function J^iiXjQ^) defined in Eq.(4.22). 

and massive resummed (solid) predictions. We observe that the two fixed-order 
calculations yield different predictions throughout the full x range, which is a 
consequence of the implementation of mass effects; however, at large x, where 
one starts to be sensitive to the resummation, the impact of soft resummation is 
competitive with that of mass contributions. 

Before closing this section, we would like to point out that, although we have 
improved our perturbative prediction by implementing soft-gluon resummation 
in the coefficient function, non-perturbative corrections are still missing. Non- 
perturbative effects are important especially at small values of and large x. 
In fact, a more accurate investigation of the very- large x limit of the structure 
functions shows that they exhibit an oscillatory behavior once x gets closer to the 
maximum value which is kinematically allowed. 
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Conclusions 

In this Thesis we have apphed perturbative QCD to make precision predictions 
for processes hke top decay and heavy quark production in charged-current DIS. 

We have first discussed the general properties of QCD as a theory of the 
strong interactions. The discussion is organized in a way that emphasizes certain 
problems that arise in processes involving heavy quarks, like the appearance of 
large quasi-coUinear logs and large soft-logs. Furthermore, in addition to the 
factorization theorem and parton evolution, we presented a general description of 
the resummation methods used to treat those problems. 

Our first application was to b quark fragmentation in top decay, since this will 
be an important observable for the study of the top quark once enough experimen- 
tal data is available. We first presented the results of our NLO QCD calculation 
for this process and we commented on its shortcomings: the need for all-order 
resummation of quasi-coUinear and soft logs. Our NLO perturbative result re- 
produces the known universal initial condition for the perturbative fragmentation 
function as well as the NLO result for the total top decay width. We performed 
the resummation of the quasi-coUinear logs ln(m^/m() applying the formalism 
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of perturbative fragmentation function with NLL accuracy. The resummation of 
those logs leads to very serious improvement of the 6-quark spectrum. To further 
improve the behavior of the 6-spectrum towards the region of large 6-momentum 
fraction, we have evaluated the soft-gluon resummed expression for the top-decay 
coefficient function with NLL accuracy. The combined result of the two types of 
all-order resummations that we have performed is significant: the validity of our 
result is extended towards higher values of the 6-quark energy fraction Xf,- The 
resulting expression has very little dependence on the unphysical factorization 
and renormalization scales. Finally, we used our improved perturbative result for 
6-quark fragmentation to make a prediction for the energy spectrum of certain 
6-fiavored hadrons in top decay. To that end we extracted information about the 
non-perturbative transition 6-quark to 6-flavored hadron from e'^e" collision data. 
We performed that extraction in two ways: first we extracted the non-perturbative 
fragmentation function by fitting e'^e~ data with the fragmentation function of 
three different functional forms having one or two fittable parameters. We have 
observed that working this way, the largest uncertainty in the prediction comes 
from the fits to the data. The second approach we pursued in order to extract in- 
formation for the non-perturbative fragmentation, is to extract from data the four 
lowest moments of the fragmentation function by working directly in the space 
of Mellin moments. That way we do not have to extract the non-perturbative 
function, i.e. no fit with its corresponding uncertainty is involved. 

The second application that we considered in this Thesis is heavy quark pro- 
duction in inclusive CC DIS. We started by making the observation that the NLO 
light quark initiated coefficient function for heavy quark production in CC DIS 
contains terms that become arbitrarily large towards high values of the Bjorkcn 
variable. That fact signals the need for soft-gluon resummation for that coefficient 
function. We carefully identified the sources of such large contributions. Then 
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we performed the soft-gluon resummation separately in the kinematical regimes 
of large values of the ratio of the mass of the heavy quark with respect to the 
transferred momentum and small values of that ratio (i.e. the massless regime). 
Our calculation reproduces the known results in the limit of zero mass and of 
the NLO evaluation of the coefficient function. We applied our result to charm 
production in CC DIS. The most relevant application is to charm production in 
neutrino-nucleon DIS at NuTeV, since that experiment can reach low values of 
momentum transfer where our resummed result is most relevant. We made a pre- 
diction for the charm-component of the structure functions with full account for 
target mass effects. We observed that at large Bjorken x, although the absolute 
value of the structure function is small, the relative impact of the resummation we 
performed is significant. The dependence on factorization/renormalization scales 
is reduced as a result of the resummation. We have checked that at low Q, the 
size of the effect of the threshold resummation (at large x) is comparable or even 
larger than the effect of the power corrections of the charm mass. This is clearly 
another indication of the significance of our result. It will be very interesting to 
compare our prediction with the forthcoming data from NuTeV. 
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Appendix A 

Some Supplementary Results 



A.l Relevant Feynman Rules 



We can specify the Feynman rules that we have used by presenting the relevant 
amplitude for one of the diagrams with real gluon emission. 



2 ' 



Cl,C2 



X 



-t^Vt,S,,,e,r (i - 7') 



{Pb+Pg^ -ml + ie 
^l{Pw,p)ux,ci{pt)- 



(A.l) 



The terms in square brackets correspond to the boson-quark vertexes, 01,2,3 are 
the quark color indexes, cr, p and A are the polarizations (spins) of the 6-quark, 
the gluon, iy+-boson and the top-quark respectively. Also, the propagator for a 
gluon with loop momentum q is (we choose to work in Feynman gauge): 



—iS, 



ab 



(f + ie 



(A.2) 



We work with cross-section that is summed over the final states and averaged over 
the initial ones. The relevant results (for the example of Eq.(A.l)) are: 



XIX] o X] X/ oX^X^X/X/' 

in out ^ ci,C2 a ^ \ a p k 
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(A.3) 
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Y^Gt{p,K)Gl{p,K) = -g,„ (A.4) 

K 

^U„,a{p)Ua,l3{p) = {^ + m)al3, (A.5) 
cr 

^liP^ p)^Ap, p) = -9v,x + (A.6) 
Our conventions for working in D space-time dimensions are: 

A-D 



e — 



2 

< = ^ 

measure for the loop integrals. (A. 7) 



(27r)^ 

Also, only the strong coupling picks up the usual mass dimension 



gs gsl^^- 



A. 2 Phase Space for Top Decay at NLO 

Wc consider the three particle top decay t bWg in D space-time dimensions, 
where each of the particles in that reaction is on-mass-shell. The masses are 
rrit, mw, iTib and respectively. For simplicity, we work in the top rest frame. One 
has: 



dr = Nd \M\'^{pt;pb,Pg,Pw) S^{pt-Pb-Pg-Pw) 

X S{pl-ml)6{pl,-mlr)S{pl)d^Pbd^Pgd^pw, {A.i 
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with: 



Nd = ^ , ,\ , ■ (A.9) 



With the help of the delta-functions, one can reduce (A. 8) to: 

dr = Nd \ W S {^{Ph-Pg) - ^{Pt-Pa) - ^{Pt-Pb) + m\ + ml - 



The matrix element |Mp is given in Eq.(3.7); we have used the relations implied 
by the delta functions to simplify it. 

Next, it will be convenient to work in spherical co-ordinates with the z-axis 

along the D — 1 dimensional vector p^. Clearly, nothing depends on its orientation; 
the only nontrivial integrations are over the energies E}, and Eg and over cos 9 = 
cos {pI.pI)- 

To carry out the trivial angular integrations, one may use the following results: 
in d-dimensional Euchdean space (see also [5]), 

d'^x = r'^-^drd^d , dnd=l[ {sm0if~^~^ d0i. (A.ll) 

1=1 

The trivial angular integrations (all for the Pg case and all but one for the pi, case) 
amount to the factor Sd ^d-, where: 



J dQd: 



^d : dnd=^d{sme)''-^de, (A.12) 

i.e., the factor ^d contains all the angular integrations but the azimuthal one. 
To evaluate those factors, it is convenient to introduce: 

TT r(7/ + 1) 

\2- 



Jo 



9d9 



2" (r(n/2 + 1))^ 
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The latter satisfy the relation J^'^) = (n-l)/nj("-2) ^j^h = yr, J^^) = 2. Then 
one can easily show that Sa = 2nn=o'^^''^ and = 2n^=o'^^"^- The relevant 
result for us is: 



S 



D-1 



%-i = 2 " 



r(3/2-e)' 



r(i-6)' 

Sd-1 ^d-1 = Stt" [l + 2e (1 - 7£ - ln(27r)) + 0{€^)] . (A.13) 

Finally, we can write the final result containing the two non-trivial integrations 
over the two energies: 



where xe = Eb/rrit, y = Eg/rrit, A = x%~b-{l-XE)'^, B = 2{1-xe){s-xe), C = 
—{xe — s)^- The squared matrix element (3.7) should be evaluated using: 

cos(p6,p^) = .^.!'^. {xEy - Xe - y + s). 
\Pg\\Pb\ 

The hmits on the two variables are (see also the notation introduced in (3.4)): 

< Xe < s, 

s — Xe 

y-<y<y+; y± = 



1-xeT \Jxe - h 

To obtain the desired differential distribution dV/dxh one must integrate over y 
and use Xf, — xe/s ets implied by (3.5). In the integration over y the so-called plus 
prescription appears. It is discussed in Section (A.4). 



A. 3 Spence Functions 

We present some relations between the Spence functions (or di-logarithms) that 
we have used. Most are taken from [105] where more details about the properties 
of those functions can be found. 
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That function is defined through the integrah 

dt 



U^{z) = - -ln(l-t), (A.15) 
Jo t 



and has a cut from 2; = 1 to oo. Some particular values are: 



U,{1) = -, U,{-1) = --. (A.16) 



The Spence function satisfies the following relations [105]: 

Li.(-i) = -Li.(-x)-iln^(x)-^, 

Li,{x) = -Li2(^)-^ln2(l-a:), 
^112(0;^) = Li2(x) + Li2(-x), 



1 — xl — yJ \1 — y J \l — x, 

- Li2(x) - Li2(?/) - ln(l - x) ln(l - y), 



— h m[x) m 



6 ^ ' Vi-yy 

Combining the above relations one can derive other relations: 

Li2(x) = -Li2(l - x) + — - ln{x) ln(l - x), 
U,{x) = Li2(^)+ln(l-x)ln(l£^)-^W(l-x)-^, 
^^^{~~^) ^ 2Li2(^Y~~^) ~3Li2(l -x) -21n(a;)ln(l -x) 



+ ln'{l-x)-^\n'{x) + ^. (A.17) 
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A. 4 Plus Prescription 

Here we discuss some properties of the so-called plus prescription [/(a;)]+ defined 
through: 

[\f{x)Ug(x)dx^ f f{x){g{x)-g{l))dx (A.18) 
Jo Jo 

where g is a sufficiently regular function, while the function / typically is non- 
integrable in the point x = 1. The plus prescriptions are strictly distributions; for 
X < 1 they can be thought as the function itself i.e.: 

[f{x)]+ = f{x) for X <1. 

Those distributions appear through the following identity: 

One can prove it by multiplying both sides with a test function, use the definition 
(A.18) and then compare the terms multiplying the different powers of e. The 
lower limit can be any number x^in '■ < Xmin < 1. In the massless case, mf, = 0, 
one needs the identity (A. 19) including the term linear in e, the latter being: 

ln(l — x) \ ln(a;) 
1 — X ) ^ 1 — X 

Using Eq.(A.18), one can easily prove the following important distributional 
identities: 

[f{x)]^ g{x) = f{x)g{x) - {g{l) £ f{y)dyj 5(1 - x), 

[fix)g{x)]^ = [f{x)]^gix)-Qjf{y)]^g{y)dyy{l-x), 
[f{x)lg{x) = [f{x)lg{l) + f{x){g{x)-g{l)). (A.20) 

All of these identities are applicable to functions /, g for which the identities 
themselves make sense. As an example, using the above relations, one can write 
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the splitting function P^'^^x) in Eq.(2.34) as: 



l-x 



((T^ - (1 + + - (.A.21) 



A. 5 Mellin- Space Results 

The MeUin transformation is defined through Eq.(2.39). That transformation 
has the factorization property (2.40). The inverse MeUin transform is given in 
Eq.(2.65). The following functions are often encountered: 

1 



/' 



dz z 



-Z) + 

In z 



dz z^-i- 

(1-z 



f 

Jo 



dz z 



N-l 



1 


'ln{l-zy 


1 r 




1-z 





/ dz z^-Hn(l-x) = 
Jo 



= -Si(N-l), 

Sl{N - 1) + S2{N - 1) 

Si{N) 



10 N ' 

where (we correct a typo in that definition in Eq.(B.2) in [51]): 

dk+i \^Y{x) 



and 



S^{N) = ^o(Ar + l)-V;o(l)=-Si(iV-l) + 



TV' 



S2{N) 



-i^i{N + 1) + zAi(l) = S2{N - 1) + 



S^{N) = -[^,{N + l)-i;2{l)]. 



(A.22) 



(A.23) 



(A.24) 



The large N behavior of some of the often encountered singular functions is [53]: 

1 



^0 



dz z 



N-l 



(l-z) 



+ 



^ -\n{N)-^E + 0{l/N), 
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Jo 



dz z 



N-l 



ln(l - z) 
1-z 



-1 + 



= ^ln2(7V)+7^1n(7V) 



(A.25) 



A. 6 A^-space Result for the Coefficient Function 



and Dl""' 



The A'"-space expression for the coefficient function (3.21) is [51]: 



27r I 



+ [l + 21n(l-w)] 
4w(l — w 



1 



-2tlji{N)-2iji{N + 2) 



+ 



l + 2w 



N{N + 1) 
2F^{l,N + l,N + 2,l-w) 



N + 1 



2Fi{l,N + 2,N + 3,l-w) 
N + 2 



+ Sl{N-l) + Sl{N+l) 



+ S2{N - 1) + S2iN + 1) + 2[1 - 2 ln(l - w)]Si{N) 

1 — IV ^zv 
+ 2 In w ln(l — w) — 2 ln(l — w) — In w 



+ 4Li2(l - - 6 - 



1 + 2w 

27r2' 



\ — w 



(A.26) 



The A'"-space expression for the quark initiated initial condition (2.58) for the 
perturbative fragmentation function is [42]: 



21 



25f(iV) + 



2 ] 1^2 iV(iV + l) 



2^2 (A^) + 2 



A^(A^+ 1) 
1 



m 

Si{N) 



-2S,{N) 



+ 2Si{N) 



(A.27) 



N(N + 1) 

One can derive those transformations with the use of the results in Section (A. 5). 
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A. 7 AT-space Expressions for the Kernels Pjy 

Applying the results in Section (A. 5) one has the following A^-space expression 
for the leading order splitting function: 

The result for the non-singlet NLO kernel is reported in [42] and reads: 

Pn = C^F [MN) + A{N)] + ^CfCaPg{N) + ^CpUfPNEiN), (A.29) 
where Cp and Ca are given in (2.10), Uf is the number of the active flavors, and 
PAN) = {2Sm-j^^j^)(2S,(N)-2m)-^^^^S,(N) 

+ 453(Ar)-35.(Ar) + 3C(2)+ ^3 +^^^3 - -, 

P (N) - ^ ^ nN' + 5N-3 

Pnf{N) - -S,{N)--S,{N)---2 9^2(;v+l)2 ' 

Pc{N) = -P,{N) + S,{N) (^-igl - 2 ^,^^^+ \^, ) + 4^,(Ar)^,(Ar) 

/13 2 \ 43 151A^^ + 263Ar3 + 97Ar2 + 3Ar + 9 

+ ^2[N}\— -————] + — + 



3 A^(Ar+l)y 24 9Ar3(7V + i; 



3 



A{N) = 2f-25,(Ar) + ^ + — -L— V252(Ar)-l^^ 2^ + ^ 



2 Ar(iv + i);v 3 iV2(iv + i)2y 

A. 8 Derivation of the Factors Ki in Eqns.(4.47) 
and (4.58) 

Here we present in some detail the derivation of the factor given in Eqns.(4.47) 
and (4.58). We introduce a functional A, such that to a function f{N): 

hm f{N) = Yl^n^N) + /o + 0(1/Ar), 

i>0 
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it assigns the constant /o, i.e.: 

A[/(iV)] = /o. 
Then the desired coefficient Ki is defined as: 

To calculate the action of A, one first has to extract the singular dependence of 
a function using the results in (A. 20). For example, making use of the relations 
(A. 25), one gets for the leading order splitting function in Eq.(A.21): 

'3 



-27. 



E 



We next calculate explicitly one of the more comphcated terms in the coefficient 
function Hf that can be found in [92]: 

^(A + ln(l-A)), A<1; 
-lE, A = 1, 



(/a(^)). 



(A.30) 



where 



1 



To demonstrate that result, one needs to separately consider the cases A < 1 and 

A = 1. The action of A in the latter case is obvious. When A < 1, fx{z) is 

integrable and we extract the z ^ 1 singular contribution using (A. 20): 

1 



[h{z)]+ = h{z) + ^ (A + ln(l - A)) S{l-z), A < 1. (A.31) 

However, one must ensure that the first term in (A.31) vanishes in the large 

N limit. One can show that this term is indeed suppressed and does not 

contribute. However, for large but fixed N, in the limit A — > 1: 

ln(l - A) 



fx{z) ^ 



N 



i.e. although the first term in (A.31) is formally 1/N suppressed, it is no longer 
small in the limit of vanishing quark mass. 



